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1.  Introduction 

Isoplanatism,  as  it  is  conventionally  used,  refers  to  the  uniformity 
of  the  quality  of  imagery  over  some  portion  of  a focal  plane.  Introduced 
into  the  literature  in  this  form,  it  serves  to  "validate"  the  concept  of  the 
modulation  transfer  function.  Its  meaning  is  physically  straightforward 
and  intuitively  quite  clear.  But  the  concept  of  isoplanatism  as  applied  to 
special  imaging-through-turbulence  concepts  is  a much  more  complex 
matter.  Here,  rather  than  defining  the  focal  plane  uniformity  of  image 
quality,  it  is  directly  tied  to  the  achievability  of  image  quality. 

Our  reference  to  special  imaging-through-turbulence  concepts  is 
meant  to  include  adaptive  optics  and  the  two  speckle  imagery  techniques, 
the  first  demonstrated  by  Labeyrie1  and  the  second  conceived  and  analyzed 
by  Knox  and  Thompson.*  (Analysis  of  the  amplitude  interferometer  concept 
of  Currie3  can  probably  be  obtained  as  a special  case  of  the  speckle  inter- 
ferometer, but  we  shall  not  explicitly  concern  ourselves  with  this  here.) 

Adaptive  optics  and  speckle  imagery  have  the  highly  significant  feature 
that  they  sense  the  received  wavefront  in  time  periods  shorter  than  that 
in  which  the  atmospheric  turbulence -induced  wavefront  distortion  can  change. 
Isoplanatism  for  these  special  imaging  techniques  refers  to  the  dependence 
on  field-angle  of  the  wavefront  distortion  (rather  than  to  the  dependence  on 
field-angle  of  the  image  quality).  With  these  special  imaging  techniques, 
if  the  wavefront  distortion  is  nearly  the  same  for  a large  enough  range  of 
field-angles,  the  "ideal"  resolution  can  be  achieved,  no  matter  how  severe 
the  actual  wavefront  distortion  per  se  is. 

We  have  elsewhere4*3  presented  an  analysis  of  the  isoplanatic  con- 
siderations in  adaptive  optics,  and  therefore  will  not  treat  that  problem 
here.  In  this  work,  we  shall  restrict  our  attention  to  isoplanatic  effects 
ii.  speckle  interferometry  (the  Labeyrie  technique)  and  in  speckle  imagery 
(the  Knox- Thompson  concept).  In  the  next  two  sections,  we  shall  examine 
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the  Labeyrie  technique  in  the  absence  of  isoplanatism  effects.  This  will 
establish  the  basis  for  our  subsequent  analysis.  The  first  of  these  two 
sections,  i.e.  , Section  2,  will  treat  a single  point  source  (a  situation  in 
which  there  is  manifestly  no  anisoplanatism  effect),  and  will  closely 
parallel  the  analysis  of  Korxf.®  The  second  of  these  two  sections,  i.e.. 
Section  3,  treats  the  same  problem,  but  for  a pair  of  point  sources.  This 
gets  us  into  a situation  where  we  could  examine  isoplanatism  effects,  but 
at  this  point  we  develop  results  assuming  that  the  point  sources  are  close 
enough  to  each  other  that  there  is  no  anisoplanatism  problem.  The  com- 
pletion of  this  analysis  prepares  all  of  the  background  that  we  need,  and 
in  Section  4 we  are  able  to  analyze  the  dual  point  source  problem  with  a 
large  enough  separation  between  the  points  that  anisoplanatism  effects 
have  to  be  considered.  This  provides  us  with  the  basic  isoplanatism  re- 
sults for  the  Labeyrie  technique.  In  the  two  sections  following  this,  we 
switch  our  attention  to  the  Knox-Thompson  concept.  In  Section  5,  we 
treat  a pair  of  point  sources  close  enough  together  that  there  is  no  aniso- 
planatism effect.  The  results  presented  here  parallel  the  basic  work  of 
Knox. 7 In  Section  6,  we  extend  this  analysis  to  allow  consideration  of  the 
case  in  which  the  two  point  sources  are  far  enough  apart  that  anisoplana- 
tism effects  have  to  be  taken  into  account.  In  this  section,  we  develop  our 
fundamental  results  for  isoplanatism  for  the  Knox-Thompson  concept. 
Section  7,  the  final  section,  reviews  all  of  the  results  and  indicates  their 
significance. 


Our  approach  throughout  this  paper  will  be  to  use  a notation  that  is 
sufficiently  general  that  it  allows  treatment  of  all  of  the  problems.  This 
will  permit  us  to  use  one  self-consistent  notation  throughout  the  paper.  It 
will,  of  course,  force  us  to  utilize  an  overly  complex  notation  in  some 
places,  but  while  this  may  appear  locally  inelegant,  it  will  make  the  totality 
of  the  paper  easier  to  read. 
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Our  analysis  will  be  directed  to  the  problem  of  optical  propagation 
down  through  the  atmosphere,  and  small  to  moderate  zenith  angles.  For 
this  case,  we  may  assume  that  the  intensity  fluctuation  effects,  as  meas- 
ured by  the  log-amplitude  variance,®  are  small  compared  to  the  phase 
fluctuation  statistics  — small  enough  to  be  of  negligible  significance. 
Accordingly,  we  shall  base  our  analysis  on  a formulation  in  which  there 
is  only  a turbulence -induced  phase  variation  at  the  entrance  aperture  of 
the  telescope j and  no  corresponding  intensity  fluctuation  effects.  In  the 
same  vein  of  reasoning,  we  shall  associate  results  obtained  in  terms  of 
the  phase  structure  function  as  though  they  applied  to  the  wave -structure 
function.®  It  can  be  shown  that  the  discrepancy  introduced  by  this  approach 
is  of  the  order  of  an  exponential  of  fee  log-amplitude  variance, so  that  it 
should  consequently  be  quite  small  for  the  problem  scenario  of  interest  to 
us  here,  i.  e.  , astronomical  viewing  geometry. 


:• 


2.  Labeyrie  Technique:  Single  Point  Source 

We  consider  a single  point  source  at  an  angular  position  in  our 
field -of -view  of  . We  assume  that  the  source  intensity,  at  wavelength 

\ , is  such  that  the  amplitude  at  our  telescope  aperture  is  .If  we 
denote  position  on  the  telescope  aperture  by  the  two-dimensional  vector  "r  , 
and  let  0 (t3  ;r)  denote  the  instantaneous  random  phase  shift  at  "r  intro- 
duced by  atmospheric  turbulence,  then  we  can  write  for  the  wave  function 
at  the  telescope  aperture  due  to  this  point  source 

U5  (?)  = Aj  exp  [ih'fj  • r + i 0(^  ;r)]  , (1) 


where 


k = 2tt/\ 


is  the  optical  wave  number. 

We  assume  that  the  telescope  has  a circular  clear  aperture  with  a 
diameter,  D . It  is  convenient  to  introduce  here  the  function 


W(?)  = 


1 . if 


0 , if 


D 

>*  D 


to  provide  a mathematical  definition  of  the  aperture.  Defining  position 
in  the  telescope  plane  in  terms  of  the  corresponding  field-angle,  , the 
wave  function  due  to  our  point  source  as  it  appears  on  the  focal  plane  can 
be  written  as 

uj  (t)  = ?[  J dr  W(r)  exp  (-ik^*-?)  U3  (?)  , 

where  si  is  a constant  of  proportionality  whose  exact  value  is  of  no 
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particular  concern  to  us  as  it  will  drop  out  of  our  final  results. 
The  focal  plane  intensity  can  be  written  as 


!(?)=*  |u,(f)|* 


(5) 


The  quantities  of  immediate  interest  to  us  are  the  fourier  (angular)  spatial 
frequency  transform  of  the  focal  plane  intensity, 


S(T)  * J exp  (-Zni?.*)  I(tf) 


(6) 


where  f is  the  spatial  frequency  (in  units  of  cycles  per  radian-field- of  - 
view),  and  a quantity,  */(!*,!*')  , which  we  call  the  bispectrum,  and  which 
is  defined  by  the  equation 


jG,?')  -<S*(?')  S(f)> 


(7) 


Here  and  throughout  this  work  we  use  the  angle  brackets  to  denote  an 
ensemble  average  over  all  possible  realizations  of  the  turbulence -induced 
wavefront  distortion  effects.  It  should  be  noted  that  in  the  case  where  the 
two  spatial  frequencies,  T and  7*  , are  equal,  the  bispectrum  reduces 
to  the  ordinary  power  spectrum.  Our  notation,  i.  e.  , ,*>(?, ?')  , allowing 
for  two  distinct  spatial  frequencies  is  more  general  than  we  shall  need  for 
analysis  of  the  Labeyrie  technique,  where  7 and  ?'  will  always  be  equal. 
However,  use  of  this  generalized  notation  will  allow  a smooth  transition 
of  the  analytic  work  into  Sections  5 and  6,  where  we  will  analyze  the  Knox- 
Thompson  concept. 


It  should  be  noted  that  we  have  dropped  the  subscript,  j , from 
the  definition  of  I(r)  , the  focal  plane  power  density.  This  will  allow 
the  same  notation  for  focal  plane  intensity,  I , and  thus  for  its  transform, 
S , and  the  bispectrum,  J , to  be  used  whether  we  are  considering  a 
single  point  source,  (j)  , or  a pair  of  point  sources,  (j  and  j')  . 
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If  we  substitute  Eq.  (1)  into  Eq.  (4),  we  get 

uj(?)=||Aj  J‘d?W(r)exp[i0(T1;-?).ik(?-^1).ri  . (8) 

Substituting  this  into  Eq.  (5),  and  rewriting  the  product  of  integrals  as  a 
double  integral,  we  get 

I(V)  = * |ll |8  A»  JJ  dr  d?'  W(?)  W(?') 

x exp  -0(YJ ;■?')] 

- ikCtt-fj)*  (?-r')]}  . (9) 

In  accordance  with  Eq.  (6),  we  can  write  for  the  spatial  frequency  fourier 


transform 


s(f)  = i 111!8  A » JJJ  d?  dr  dr'  W(r)  W(?') 

X exp  (-2tt  i T*?)  exp  {i  [0(^  ;"?)  -0(tf  ;r')] 
- (r-7#)]} 


Replacing  the  variable  of  integration  J by  ? = * we  8et  *n  place  °f 

Eq.  (10) 


S(F ) = £ |u|s  A^  exp  (-ZTTif.^)  [JJ  d?  <f?  dr ' 
X W(?)  W( ?')  exp  { i [ 0 ;7)  -0(1/?')] 

- ik?.  (?+  Xr  -?')] 
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It  is  a well-known  property  of  the  fourier  transformation  that  the 
repeated  fourier  integral  recovers  the  starting  function.  Thus  for  h(y) 
being  any  reasonably  well  behaved  function 

JJ  dx  dy*  exp  [ ±Zr\  i x(y-y')]  h(y))  = h(y)  . (12) 

Making  use  of  this  result,  we  can  perform  the  and  “r '-integrations  in 
Eq.  (11).  Thus  we  obtain  the  result  that 

S(T)  =$  |«|»  A*  exp  (-2TTif.  J*  d?  W(?)  W(r+\?) 

X exp  (i[0(^  ;r) fr+X?)]}  , (13) 

where  we  have  suppressed  a factor  of  X8  in  this  result  by  lumping  it  in 
with  the  constant  | f|  | a . 

We  calculate  the  bispectrum,  or  rather  since  we  are  interested  in 
the  case  where  T = T*  , the  power  spectrum  ^(f\f)  , which  we  will 
write  as 

/(?)=^(r;r)  , (i4> 

1 1 

by  substituting  Eq.  (13)  into  Eq.  (7),  This  gives  us  the  result  that 

J CD  = \ IsiIMA,8)8  J*f  d?  d?'  W(?)  W(?+xf) 

x W(r')  W(?'+X?)  (exp  {i[0(?3;r  ) -*(?,;?') 

+ 0(t3;?'+X?)  -0(?1;?+Xf')]}>  , (15) 

where,  in  order  to  obtain  this  result,  we  have  made  a double  integral  of 
the  product  of  integrals. 


To  reduca  thi»  expression,  we  start  byr  noting  that  the  quantity  in 
the  equare  brackets  in  the  exponential  is  a gaussian  random  variable  with 
aero  mean,  and  that  for  such  a random  variable,  say  ^ , with  g being 
any  constant, 

<«xp  (cxx»  = exp  (*  o»  <xa»  . (16) 

Making  use  of  this  result  and  the  fact  that 

(a-b+c-d)*  = (a-b)»-  (a-c)»  + (a-d)»  + (b-c)«-  (b-d)»  + (c-d)« , (17) 
we  can  rewrite  Eq.  (15)  as 

~ k |»|4(A  *)»  JJ  dr  €t*  W(r)  W (r+\T) 

X W( f')  W(r'+\T)  exp  {-$  [j<?-?') 

- -\?)+Jt(\T)+£(\7)  -J}(r  -f'+\T) 

+ ^(? -■?*)]}  , (18) 

where  .0(p)  is  the  wave -structure  function,  defined  by  the  equation 

*(?>  = <C0(T,  ;"?+•?)  -*(*,;*>]•>  . (19) 

Consideration  of  the  statistical  symmetry  of  the  wave -structure  function 
makes  it  clear  that  J>(p)  = , a fact  which  we  have  used  in  developing 

Eq.  (18)  in  the  form  presented. 

It  is  obvious  that  when  we  consider  the  zero  frequency  case,  i.e.  , 
when  = 0 , the  exponential  in  Eq.  (18)  goes  to  exp  (0)  = 1 , so  that  we 


can  write 


•'’(0)  |«|*(AJ»)*  J'fdr  dr'  [W(r)]»[W(r')]» 
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Consideration  of  Eq.  (3)  makes  it  apparent  that  W-squared  is  equal  to 
W . Moreover,  we  note  that  the  double  integral  of  Eq.  (20)  can  be  re- 
written as  the  product  of  two  identical  integrals,  each  equal  to  D®  . 

Thus  we  can  rewrite  Eq.  (20)  as 

JiO)  =(£  |tr|»  A»  D»)'  . (21) 

This  will  provide  a basis  for  normalization  of  the  results  we  shall  develop 
shortly. 

To  proceed  with  the  reduction  of  Eq.  (18),  we  make  use  of  the  fact 
that  the  wave  - structure  function  can  be  written  asii.18 

2>(p)  = 6.88  (p/r0)Bfe  . (22) 

where  rQ  is  a length  which  serves  as  a measure  of  the  transverse  co- 
herence of  the  turbulence  distorted  wavefront.  Thus  we  can  rewrite  the 
exponent  in  Eq.  (18)  as 

E = Ofr-r ')  - + j(\?)  + j>(X?) 

- Jfr(r-r'+  \T)  + jtfr-r ')] 

= - [M r-r')  + M\T)  - i X?)  - £ M ?-?'+  Xf )] 

= - 6.  88  rQ-6/&  ((?-?' |e/&  + | I?--?'-  X?|e* 

- £ |r-r'+Xf  |6/&)  . (23) 

Based  on  the  approximation  that 

l + f6/s.|(l  + 2e  cos  9 + c3)*'8  (1  - 2c  cos  e + e3)B/6 

« [i-I  ^(1-Jcob*  8)]  . (24) 
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which  is  accurate  for  small  enough  values  of  « , and  noting  that,  by  virtue 
of  die  law  of  cosines,  we  can  write 


t 

I 


|r-?'±\f  |s*  = (|r--?'|*  ± 2|"?-‘?/|  |\? | cos  . (25) 

where  9 is  the  angle  between  the  vectors  (r-"r')  and  (x7)  , we  can  in- 
troduce die  approximations  that 

E«-6.88  (X|T|/r0)e/3|-l.|  (-^jpt-^cos*  9)]  . 

if  X|?|  « jr-"r'|  » 

and 

E *«  - 6.  88  (|^-?^|/r0)6^|-l-|  cos*  9)]  , 

if  « ,X|^|  • (27) 

In  developing  Eq.  (26)  from  Eq.  (23),  terms  of  higher  order  in  | X? | / 1"?— "?  ' ] 
have  been  dropped  from  inside  the  square  brackets.  [The  first  term  dropped 
is  proportional  to  (X)?|  / |?-r ' | Similarly,  in  developing  Eq.  (27)  from 

Eq.  (23),  terms  of  higher  order  in  |‘r-r'|/|x7|  have  been  dropped  from  in- 
side the  square  brackets  on  the  right  hand  side  of  Eq.  (27).  [Here  the  first 
term  dropped  is  proportional  to  (I”?--?'!  /x|7|  )7/3.  ] For  our  purposes,  it 
will  be  sufficient  to  replace  Eq.  's  (26)  and  (27)  with  the  further  approxima- 
tion that 


E~-  6.88  (x|7|/r0)fi*  , if  X|fj  « \*-?'\  , (28) 

E w-6. 88  ()?-?' | /r0)e/&  , if  |r--?'|  « X|f|  . (29) 

If  we  study  the  exponential  in  Eq.  (18)  subject  to  the  assumption 
that  X|f|  is  very  large  compared  to  rQ  , we  see  that  the  only  condition 
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under  which  the  exponential  function  will  not  be  vanishingly  small  is  when 
| | is  not  much  larger  than  (or  smaller  than)  r0  . In  this  case,  in 
accordance  with  Eq.  (29),  we  can  write 

J(f)*$  l«!4  (A  »)•  J*J*  dr  dr'  W(r)  W(r+\T)  W(r')  W (r'+X?) 

X exp  [-6.88  ( | "r "|/ r0 je/a ] . (30) 

It  will  be  noted  that  the  double  integral  in  Eq.  (30)  appears  to  allow  the 
case  where  |‘?-'r/j  is  not  particularly  small.  Strictly  speaking,  we  should 
include  some  condition  in  the  integrand  restricting  the  range  of  integration 
so  that  we  only  considered  pairs  of  values  of  (r, ~r*)  for  which  |’r-r'|  was 
comparable  to  or  smaller  than  rQ  . In  fact,  however,  such  a restriction 
is  automatically  provided  by  the  exponential  function  in  Eq.  (30). 

Based  on  the  assumption  that  the  aperture  diameter  D is  much 
larger  than  rfl  , and  noting  that  the  form  of  the  exponential  makes  the 
integrand  vanish  when  ~r  is  not  nearly  equal  to  ~r*  (i.  e.  , when  j’r-'r'  | 
is  not  comparable  to  or  smaller  than  rQ  ),  we  can  make  the  approximations 

that 


W (r)  « W(?') 

in  Eq.  (30)  , 

(31) 

Wf?+X?)  * W(?'+Xf)  , 

in  Eq.  (30) 

(32) 

In  accordance  with  Eq.  (3),  we 

can  write 

W(r)  W(?')  ® W (?) 

, in  Eq.  (30)  . 

(33) 

W(r+Xf)  W(r'+X?)  « W (r+X?)  , in  Eq.  (30)  . (34) 

If  we  make  a change  of  variable  of  integration  in  Eq.  (30)  from  (r,  r') 
to  (r,x)  , where 
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X = r-7' 


(35) 


then  with  the  use  of  Eq.  's  (33)  and  (34),  we  can  rewrite  Eq.  (30)  as 

S<T)  = * fill*  (Aj®)*  { j d?  WC?)  W(?+Xr)| 

X j J*  dx  exp  [-6.  88  (|x| /r0)B/h]|  . (36) 

The  integral  in  the  second  curly  brackets  can  be  evaluated  by  rewriting  it 
in  polar  coordinates,  performing  the  azimuthal  integration,  and  making  a 
change  of  variable  for  the  radial  integration  so  as  to  cast  the  integral  in  the 
form  of  a gamma -function,  T(6/5)  to  be  exact.  Thus  it  can  be  shown  that 

J*  dx  exp  [-6.  88  (|x|  /r0)s/3 
00 

= 2tt  J x dx  exp  [-6.  88  (x/^)®/®] 
o 

CO 

= [tt/(6.  88/r0)s/6  ] J du  exp  (-u®/®) 

0 

00 

= | [tt/ (6.  88/rj®/3)0/®  ] J v*/®  dv  exp  (-v) 

o 

= | [>/(6.88/r06/3)8/6  ] r(6/ 5) 

= 0.  342  r0®  . (37) 

The  integral  in  the  first  curly  brackets  in  Eq.  (36)  can  be  recognized  as 
corresponding  to  the  area  of  overlap  of  two  circles  of  diameter  D , with 
center-to-center  separation  of  X?  . This  corresponds  to  the  aperture 
area,  1 tt  D®  , times  the  diffraction-limited  optical  transfer  function, 
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(38) 


T0l  (?)  = (coffi(\|r|/D)  - (X|T|  /D)  [1  -(X|T|/D)8]l/S} 


for  optica  of  diameter  D operating  at  wavelength  X . Thus  Eq.  (36)  can 
be  rewritten  as 


.?(?)  = *!tl|MV)*  {*tt  EP  T0L(f)}  {0.  342  rQ*} 


With  a bit  of  manipulation  of  terms,  and  making  use  of  Eq.  (21),  this  result 

can  be  rewritten  in  normalized  form,  representing  the  speckle  transfer 
function  as 


- 0-«5  (r0/DC  Ttl(T) 


for  X|?|  » rQ 


(It  can  be  shown  that  0.43  5 is  actually  an  approximation  for  2-s/s  *«  0.  435.  ) 

If  we  now  study  Eq.  (18)  subject  to  the  assumption  that  \|f  | is 
small  compared  to  rQ  , then  we  can  see  that  in  this  case  the  exponential 
function  will  only  deviate  from  unity  when  |'r-?'|  is  large  enough  that  the 
small  difference  between  ~r-~rf  and  ~r-r'±\f  makes  the  cancellation  of 
the  £(?-■?')  and  ^(f-?'±\f)-terms  significantly  inexact.  In  this  case,  we 
may  use  Eq.  (28)  to  allow  us  to  rewrite  Eq.  (18)  as 

J(T)  |«|4(AJ*)J  JJ*d?d?'  W(?)  W(?+\?)  W(?')  W(?'+Xf ) 

X exp  [-6.88  (x|?|/r0)B/3]  . (41] 

In  this  case,  the  double  integral  can  be  decomposed  into  the  produce  of 
two  identical  integrals,  and  we  can  write 

j(T)  = {fc  |?l| * (A»)»  exp  [-3.44  (X|f  |/r0)®*] 

X J dr  W(r)  W(r+X?)^*  , (42] 
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where,  since  the  exponential  is  not  a function  of  the  variable  of  integration, 
we  have  taken  it  out  of  the  integral  (and  then  taken  its  square  root,  along 
with  that  of  £ |»|4(AJ*)S  , thus  replacing  6.88  by  3.44  ).  The  integral 
can  be  recognized  as  the  area  of  overlap  of  two  circles  of  diameter  D with 
center-to-center  separation  of  \T  . As  noted  before,  this  integral  can  be 
written  as 

f dr  W(r)  W(?+X?)  = k " I*  t01(?)  . (43) 

where  -rDL  , as  defined  by  Eq.  (38),  is  the  optical  ti-ansfer  function  for  a 
circular  aperture  optical  system  operating  at  wavelength  X . Making  use 
of  Eq.  's  (21)  and  (43),  we  can  rewrite  Eq.  (42),  in  terms  of  the  speckle 
transfer  function  as 

-JW  = [V^exp  [-6.88(x|r|/ro)B/33,for  X|f|  « rQ . (44) 

Eq.  's  (40)  and  (44),  with  Eq.  (38)  supporting,  represent  our  basic 
results  for  the  analysis  of  the  performance  of  the  Labeyrie  technique  opera- 
ting against  a single  point  source  — though  the  fact  that  these  results  are 
applicable  to  the  Labeyrie  technique  is  not  immediately  apparent  from  what 
we  have  said  thus  far.  We  shall  discuss  this  applicability  shortly,  but 
before  we  do  this,  it  is  worthwhile  to  consider  the  accuracy  of  the  approxi- 
mations we  have  made  in  developing  Eq.  's  (40)  and  (44)  from  Eq.  (18) 

[and  Eq.  (21)].  This  will  be  of  value  in  our  later  analysis  where  we  shall 
want  to  make  the  same  kind  of  approximations. 

To  demonstrate  the  accuracy  of  Eq.  's  (40)  and  (44),  we  have  carried 
out  numerical  evaluation  of  the  double  integral  in  Eq.  (18),  and  obtained 
exact  results  for  the  speckle  transfer  function,  0)  . The  integral 

evaluations  were  carried  out  using  Monte  Carlo  techniques.*  In  Fig.  's  1 to  5, 

* To  insure  reasonably  rapid  convergences  of  our  Monte  Carlo  results,  we 
used  uniform  sampling  over  the  ? and?'  circles  when  xl?|  is  less  than 
*0  — but  when  >T  is  greater  than  rfl  , we  made  the  variables  of  integration 

(Continued  on  next  page) 
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we  show  these  results  for  D/r0  * 

Also  shown  are  the  approximate  results  calculated  from  Eq.  's  (40)  and 
(44).  As  can  be  seen,  these  approximate  results  are  quite  accurate  in 
their  respective  limits — and  accordingly,  we  argue  that  the  approximation 
technique  we  used  to  develop  these  equations  is  basically  sound.  In  parti- 
cular, we  argue  that  our  method  of  inspection  of  the  exponential  in  Eq.  (18) 
to  decide  what  relationship  between  ~r-~r'  and  \f*  dominated  the  result, 
and  thus  to  justify  use  of  either  Eq.  (28)  or  Eq.  (29),  leads  to  correct  re- 
sults in  the  limits  of  }X  much  larger  or  much  smaller  than  rQ  , and  we 
shall  use  it  in  the  following  analysis. 

At  this  point,  we  are  ready  to  take  up  the  question  of  the  relation- 
ship between  the  quantity  we  have  been  evaluating,  i.  e.  , jf  (f  ) /,/  (0)  . and 
the  performance  of  the  Labeyrie  technique  for  speckle  interferometry.  In 
the  Labeyrie  technique,  the  image  of  an  object,  I(?)  , is  formed.  This  is, 
of  course,  a randomly  distorted  short  exposure  image.  The  fourier  trans- 
form of  this  image,  S(T)  is  calculated,  and  is  also  a random  function,  as  is 
I(?)  . A series  of  these  are  generated  and  represent  a randomly  sampled  set 
from  the  ensemble  of  all  possible  atmospheric  turbulence -induced  distor- 
tions of  the  image.  Using  this  set  of  image  samples  and  the  fourier  trans- 
forms, the  Labeyrie  technique  then  forms  the  (ensemble)  average  of  the 
square  of  the  transform,  ^(T.f  ) = j (t ) . This  quantity  is  taken  to  repre- 
sent a measure  of  the  source  object's  power  spectrum.  To  the  extent  that 
high  spatial  frequencies  are  present  in  /(f)  , the  correlation  function  which 
we  can  obtain  from  this  power  spectrum  will  contain  information  on  the  fine 
details  of  the  object.  In  this  sense,  Eq.  (40)  is  critically  important  since  it 
means  that  image  details  will  be  present  at  the  diffraction  limit.  The  strength 

(Continued  from  previous  page)  "r  and  . In  this  case,  we  used  uniform 

sampling  of  "?  over  its  circle,  but  sampled  ~r-~r*  according  to  the  two- 
dimensional  gaussian  distribution  with  probability  density  [2tt(0.207  r0s)]-t 
X exp  {-£  [ I ■?-■?'  | /(0.  207  r0)]3}  over  the  infinite  (?-"?' )-plane.  [We  note 
that  £ (0.  207)-e/3  *<6.88  . ] 


of  these  details  is  attenuated  by  a factor  of  (r0/D)[T6l(f  )]~1/8  compared 
to  what  would  be  obtained  with  diffraction-limited  images,  i.  e.  , in  the 
absence  of  turbulence  effects.  But  this  attenuation  factor  is  not  so  large 
as  to  be  intolerable. 

Our  analysis  to  this  point,  paralleling  that  of  Korff®  , demonstrates 
that  with  the  Labeyrie  technique,  the  high  spatial  frequencies  of  the  image 
are  not  attenuated  by  more  than  a factor  of  (r0/D)[TOL(f  )]~l/a  , and  is  of 

great  general  interest  as  it  shows  that  the  high  spatial  frequency  details 
are  recoverable.  However,  since  our  analysis  deals  with  only  a single 
point  source,  it  is  unable  to  accommodate  any  isoplanatism  effects.  To 
incorporate  isoplanatism  effects  in  our  analysis,  we  need  to  consider  a 
source  with  some  spatial  extent.  The  simplest  such  source  is  a pair  of 
point  sources.  In  the  next  section,  we  shall  present  an  analysis  parallel- 
ing the  one  in  this  section,  but  for  a pair  of  point  sources — with  the  assump- 
tion that  there  is  no  anisoplanatism  effect.  This  will  then  serve  as  a ref- 
erence for  the  analysis  in  Section  4 of  the  Labeyrie  technique  with  a pair 
of  point  sources  and  anisoplanatism  effects. 
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3.  Labeyrie  Technique;  Two  Point  Sources 

We  shall  consider  here  a pair  of  point  sources  with  angular  positions 
and  ^ in  our  field-of-view,  and  amplitudes  A and  A at  the  tele- 
scope  aperture.  We  associate  atmospheric  turbulence -induced  wavefront 
distortions,  0 ;"r)  and  0(0  ;"?)  with  the  wavefront  sensed  at  "r  due  to 

each  point  source.  Our  assumption  that  in  this  case  there  is  no  anisoplana- 
tism  problem  can  be  translated  into  the  requirement  that  0 (8  ^ ;"?)  and  0(9^',?) 
are  virtually  equal,  but  in  anticipation  of  the  problem  to  be  treated  in  Section  4, 
where  there  is  an  anisoplanatism  problem,  we  shall  retain  the  distinct  nota- 
tion for  0 (e 3 ; r ) and  0 Qf yir)  . 

We  may  write  for  the  wave -function  at  the  optics  aperture  due  to  these 
two  point  sources,  in  correspondence  with  Eq.  (1), 


Uj  (?)  © Ur(r)  = exp  [ik^  • r + i 0 (f } ;?)] 

© exp  [ikfy  + i0(fr;r)]  . (45 

The  special  notation,  © , is  used  here  to  indicate  a summation,  but  of  an 
unusual  type.  Its  proper  understanding  requires  that  we  recall  that  the 
wave -function  is  actually  subject  to  a temporal  oscillation  at  a frequency  kc 
(of  the  order  of  4x101*  rad/sec),  and  that  Uj  (?)  is  just  a phasor  assoc- 
iated with  this  oscillation,  which  we  use  to  take  note  of  the  relationship  of 
the  phase  of  the  oscillation  at  two  positions,  ? and  ~r * , by  means  of 

the  product  U*(r)  U (?')  . Moreover,  the  oscillation  itself  is  "somewhat 
erratic,"  representing  black  body  radiation  over  some  non-negligible  spec- 
tral bandwidth.  Despite  this  erratic  nature,  the  phasor  (?)  nonetheless 
represents  the  relationship  between  the  phase  at  "r  and  ?'  . However, 
there  is  no  correlation  between  the  erratic  oscillations  of  the  radiation  from 
the  two  point  sources,  denoted  by  j and  j*  , which  are  statistically  inde- 
pendent black  body  radiation  sources.  Hence  there  is  no  well-defined 
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phase  relationship  between  the  two  wave -functions,  and  U (r)  and  U^fr) 
have  no  special  meaning  with  respect  to  each  other  — their  product  has  no 
particular  meaning,  and  will  vanish  on  the  average.  Thus  we  use  our  nota- 
tion + to  denote  the  presence  of  two  statistically  independent  phasors,  with 
the  understandingthac  the  average  of  the  product  of  the  two  phasors  will 
vanish. 

In  parallel  with  Eq.  (4),  we  can  write  for  the  wave -function  at  the 
focal  plane  of  the  optical  system 

Uj  (f ) u^fif)  = 91  J*  d"?  W(r)  exp  (-i  kHf  • r ) [U^  (r ) 0 U^(r)]  . (46) 

In  parallel  with  Eq.  (5),  we  write  for  the  focal  plane  intensity 

I(fr)  = £ |Uj(J)  ft)  uJ#(J)|»  . (47) 

In  practice,  the  time  required  to  measure  the  intensity  is  so  large  com- 
pared with  the  inverse  of  the  spectral  bandwidth  of  the  black  body  radiation 
(even  for  a very  short  exposure)  that  measurement  of  the  intensity  provides 
sufficient  averaging  that  the  cross-product  u *u  t in  Eq.  (47)  will  vanish. 
Thus  we  can  write  in  place  of  Eq.  (47) 


| 


i 


I(?)  = | [!Uj  (f)|8  + !ur(?)|3] 


(48) 


With  I ($)  for  a pair  of  point  sources  written  in  this  form,  Eq. ' s (6)  and 
(7)  remain  applicable  here. 

We  can  now  write  in  place  of  Eq.  (9) 

I(?)  = * !«|a  fj*  d?  d?'  W(?)  W(r ') 

x |A3S  exp  - ik  [(?-?,)•(?-?')]  } 


t 


18 


Substituting  this  into  Eq.  (6)  and  interchanging  order  of  integration  as  appro- 
priate, we  get  for  the  random  amplitude  at  spatial  frequency  7 , 

S(?)  = £ <*9  dr*  d?'  W(?)  W(?')  exp  (-2*i?.?) 

* fAj3  exp  - ik[C?-^j) • (?-■?')]} 

+ Aj»  exp  riC0^;?)-0(fr;?')]-ik[?-^)*(?-r')]}|  . (50) 

By  separating  this  expression  into  a sum  of  integrals  and  then  replacing 
the  variable  of  integration,  , by  in  one  of  the  integrals,  and 

by  *9  = in  the  other  integral,  and  then  combining  the  sum  of  integrals 

into  a single  integral  on  a sum,  we  can  obtain  in  place  of  Eq.  (50)  the  re- 
sult that 


S(f)*i  1 91 1 3 Jf  f d?  d?  d?'  W(?)  W(r ') 


x fA38  exp  ^~2n  exP  (i  [0(?j ;“?) -0&J ;"?')]  - ik£.  (r+\f -■?')} 

+ Aj*  exp  (-2tt  if-?r)  exp  (i  [0  @y;r)  - 00^,  ;r  ')]  - i k£  • Cr+\7-'r  ')}]i 

(51) 

Now,  if  we  make  use  of  Eq.  (12),  we  can  perform  the  and  ? "-integrations. 
Thus  we  can  write 


I 

I 


S(f)  = g 1 91 1 3 J d?  W(?)  W(?+\f  ) 


x f AJ3  exp  (_2tt  ^j)  exP  {i  [0  flfj  ;■?) -0(9  J ;"r+\f)] } 

+ A»  exp  (-2tt  if*^J  exp  {i[0(er;r)-0(e’j/;r+K?)]}|.  . 
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(52) 


nMMMtt 


We  recall  here  that  just  as  0 is  a turbulence-driven  random  function,  so 
is  S(T)  a random  fourier  transform  amplitude. 


} 

1 1 

The  corresponding  power  spectral  density  associated  with  the  spatial 
frequency  t can  be  obtained  by  substituting  Eq.  (52)  into  Eq.  (7).  Making 
use  of  the  notation  simplification  introduced  by  Eq.  (14),  making  a double 

at  ] 

integral  of  the  product  of  integrals,  and  interchanging  the  order  of  ensemble 
averaging  and  integration,  we  write 

• i 

J(t)  = i | SI | 4 J'J’  dr  dr  * W(?)  W{r+Xf*)  W(r')  W(r'+\T) 

X ( (Aj*  )3 e xp  fi[0(^  ;"?) -0(f3  jr+X?) +0  (^  yr'+yj) -00^  I-?')]} 

+ (A9)(A^)  exp  [-2rri T-  (9} 

X exp  { i [0 ;r) -0^  ;r+x?)+0(fr;r#+x?)-0(^r;r')]} 

+ (A»)(AyP  exp  [2n  it . (^  -?,,)] 

X exp  {-i[0(S>J;r)  - 0 } ;(t+\T ) + 0^/?  '+\T ) - 0 ;r‘ ')] } 

+ (A  *?>  exp  {i[0(f^;?) -0(^J>;r+x?)+0(e^;?#+X?) -0(^j,;r')]}>.  (53) 

It  is  to  be  noted  that  at  this  point  we  have  retained  terms  corresponding  to 
"interference"  between  the  two  point  sources  — but  in  this  case,  it  is  "inter- 
ference"  of  the  two  random  intensity  patterns. 

The  ensemble  average  of  the  sum  of  four  terms  in  Eq.  (53)  can  be 
replaced  by  the  sum  of  four  ensemble  averages.  Following  the  procedure 
used  to  get  from  Eq.  (15)toEq.  (18),  i.  e.  , making  use  of  Eq. 's  (16)  and 
(17),  and  the  definition  in  Eq.  (19),  we  can  show  that  the  first  and  last  of 
the  ensemble  averages  can  be  reduced  to  the  form 
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(54) 


< exp  {iC0^;r)  -0(eJ;r+xf’)  + 0(^3  ;?'+*?) -0Cf,  ;■?')]}> 

= exp  ^-{^(r-r')+^(\7)  -$  CMr-'r'-Xl*) +.#(?-'? '+X?)] } ^ . 

for  the  first,  and  similarly  for  the  fourth 

( exp  {i[0(?3,;?)  -0(?J,;7+xf>)  + 0(STr;"r'+xf ) -00^,;"?')]  }> 


= exp  (-{A(r-r')  + />(X?) - £ |>(r-r '-XF)  + ^(?-r'+\D]}>  . (55) 

Here  the  results  are  essentially  identical,  since,  in  accordance  with  Eq.  (19), 
the  wave  - structure  function  does  not  depend  on  the  distinction  between  0^ 
and  fj,  . 

For  the  reduction  of  the  second  and  third  ensemble  averages  ob- 
tainable from  Eq.  (53),  the  procedure  is  more  complex.  For  example, 
for  the  second  of  the  ensemble  averages,  in  accordance  with  Eq.  (16), 
we  would  be  concerned  with 

< {[0CT3 ;?)  -0(f3  ?+XT)]  - [0  tf3,;-?')  -0(*/?'+xr)]}»  > 

= <[*(*, ;?)  -0(T3;r+xr)]3  > + <[0(9^;?')  -0(*3,;7'+xr)]»  > 

- 2 <[0(e3  ;r)-0(?3  ;7+x?)][0(03,;r')-0(r3,;?'+x?)]> 

= 2^(Xr)-2<[0^;?)-0(?i;r+\f)][0(?j,;r')-0(fr;r'+Xr)]>  . (56) 

Here  we  have  made  use  of  Eq.  (19),  but  replaced  -\T  by  \T  , based  on 
the  fact  implicit  in  Eq.  (22)  that  the  wave -structure  function  is  only  depen- 
dent on  the  magnitude  of  its  vector  argument.  With  a bit  more  algebraic 
manipulation,  Eq.  (56)  can  be  rewritten  as 
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< {[0OT, ;r) -<t>6s  $+\T)l  - ')  -t(9y;?'+\f  )]}•  > 

= z M\f  ) + ([</> 6 1 if)  -0(?J>;?')]8  > - < C0  (e , ;f)  -0(0,,  ;?'+>/)]•  \ 

-<[0(eJ;r+xf')-0(ej>;r')]a>  + <[0(fl,3;r+x?)-0(?3,;?'+xf)]a>  . (57) 

At  this  point,  making  use  of  Eq.  (16)  to  provide  a basis  for  our  use  of 
Eq.  (57)  in  the  evaluation  of  Eq.  (53),  and  substituting  Eq.  's  (54)  and 
(55)  into  Eq.  (53),  we  get 

J(t)  = * |a|*  J7  dr  dr'  W(?)  W (r+\?)  W(?')  W(?'+\f) 

X |[(A3»)»+(A3S)a]  exp  (-f^(?-r')+^(xn-^[^(?-r'-x?) 

+ ^(r-r'+x?)]}) 

+ 2(A3a)(A3*)  cos  [2n?- (03 -f3J]  exp  (-^(xf) 

' i < [0  (®  J ) - 0 (e  3/,r  ' )]®  > - ^ < [0  (f  5 ; r+  \7)  - 0 (e y fr '+x? )]a  ) 

+ ^<[0(93;?)-0(^r;r'+X^)]a>  + |<[0(93;r+xf>)-0(f3,;r')]a>}|  . (58) 

We  first  of  all  note  here  that  by  virtue  of  the  stationarity  of  the  wavefront 
distortion  statistics 

([0  (£j  ;r) -0(e3,;r')]a  > = <[0(f3  ;?+xO -0(f3,;r '+x?)]a  > . (59) 

At  this  point,  we  introduce  the  approximation  that  we  can  replace 
<[0(e  ;r)  - 0(e'3,;r')]a>  by  <0(e3;r')  - 0(93,;r)3a>  . This  approximation 

should  be  quite  accurate  when  -03,  is  very  small  (as  it  always  is  in  cases 
of  interest  to  us,  and  particularly  accurate  within  the  integral  where  all 
values  of  (r,r')  are  considered.  With  this  approximation,  and  making 
use  of  Eq.  (59),  we  can  rewrite  Eq.  (58)  as 


: 


j(l)  = i | «|*  JJ  dr  d?'  W(r)  W(?+xf)  W(r')  W (r'+\T) 

* |[(A*)«  + (Aj»)*]  exp  (-f^(?-?')+^(x7)-^[^(r^'-\r) 

+ £(?-■? '+xf)]}>  + 2(A1*)(AJ»)  cos  [>f  • dfs  -^)]  exp  {-^(xf) 

- h < [0  3 : ? ) - 0 (3 ,, * ' )38  > - * < [0  , ;?' ) - 0 (?,,  ; ?) ? > 

+ k ( [0  (fj  ;*>  - 0 0f3, ;?'+xf )]a  > + i < [0  (*3  ) - 0 (? 3, ;r +xf )]■  > 

+ i<[0(?j;r+xf)-0arj,;?')]>>  + i<[0(?j  ;?'+xf)-0(rj,;?)]8>}J 

(60) 

We  have  been  careful  up  to  this  point  to  retain  the  distinction  between 
0(?jjx)  and  0(?3,;x)  in  our  notation.  This  will  make  Eq.  (60)  a valid  start- 
ing point  for  our  analysis  in  the  next  section,  when  we  wish  to  consider  the 
effects  of  anisoplanatism.  But  in  this  section,  we  are  restricting  our  atten- 
tion to  the  case  where  {f3  and  are  sufficiently  close  together  that  we 
can  assume  that  isoplanatism  conditions  are  satisfied.  In  this  case,  we  can 
write  for  the  wave- structure  function,  in  place  of  Eq.  (19), 

~ < [0  (f3 ; r+7)  - 0(er;r)]a)  . (61) 

Making  use  of  Eq.  (61)  and  of  the  fact  that,  in  accordance  with  Eq.  (22), 

^(p)  = M-~p)  , (62) 


we  can  rewrite  Eq.  (60)  as 

«/(?)  = J 1 3I|  4 fj  dr  dr'  W(?)  W(?+xf)  W(?')  W(r'+xf) 
x |[(A»)a  + (Aja)a]  exp  ( - {*(?-?')  + ) - * [ 

+ fi(r -?'+Xf)]})+  2(A3a)a  + (A3»)a  cos  [2rt  f • (e 3 -f3>)] 


7 
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* exp  (-{^(\?)  + £(r-r')-$[^(r-r'-\?)+.0(r-r'  nr>]))J 

= i |9»|4{(a38)8+2(V)(aj*)8  cos  rzTrf.^-^jD  + tv)*} 

X JJ  dr  dr'  W(r)  W(?+Xf*)  W(r')  Wfr'+xf) 

X exp  (-U(?-?')+^(Xf)-^0(7-7'-X?)+^(?.?'+x?)]})  . (63) 


At  this  point,  we  note  that  this  expression  is  identical  to  the  one  we  had  to 
evaluate  in  the  preceding  section  for  the  single  point  source  case,  as  can 
be  seen  by  comparing  Eq.  's  (63)  and  (18).  Only  the  coefficient  (A  a)a  in 

3 

Eq.  (18)  has  been  replaced  by  the  more  complicated  expression  {(A^8)8 
+ 2(A  8)(A  8)  cos  [2nf  • ($  -0  )1  + (A  a)®  1 « but  this  is  simply  the  difference 

in  the  spatial  frequency  power  spectra  of  the  single  point  source  and  of  the 
pair  of  point  sources.  The  effect  of  atmospheric  turbulence  on  the  ability 
of  the  Labeyrie  technique  to  obtain  information  on  the  high  spatial  frequency 
part  of  the  image  power  spectrum  is  the  same,  whether  the  object  is  a 
single  point  source  or  a pair  of  point  sources,  provided  that  the  point 
sources  are  close  enough  together  that  isoplanatism  applies.  For  the  low 
spatial  frequencies,  for  an  isoplanatic  pair  of  point  sources,  the  normal- 
ized image  spatial  frequency  power  spectrum  obtained  by  the  Labeyrie  tech- 
nique can  be  inferred  from  Eq.  (44)  to  have  the  form 


JL ILL  _ 
^(0) 


(Aj8)2+  2(Aj8)(Aj8  ) cos  [2n? • -*3J]+  (A38)a 

[(a38)  + (a38)]8  KL<f>}9 


for  X | f | « r 


For  high  spatial  frequencies,  we  can  see  from  consideration  of  Eq.  (4  0) 
that  the  normalized  image  spatial  frequency  power  spectrum  obtained  by 
the  Labeyrie  technique  has  the  form 
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1 


.ZILL  « 
J(0) 


(A;)*+  2(A,»)(AJ»)  cos  [Ztt 7 • af  -T  )]  + (A,*F  r ■ . 

[(A»)  + (As*  )]»  {°-435(d)  T*L(f)} 


for  Xlf  I » r 


With  these  results  in  hand,  and  Eq.  (60)  available  for  future  use,  we 
are  now  ready  to  turn  our  attention  to  the  evaluation  of  atmospheric  turbu- 
lence effects  on  the  performance  of  the  Labeyrie  technique  in  the  absence  of 
isoplanatism.  In  the  next  section,  we  take  this  up  for  the  case  of  a pair  of 
point  sources. 
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Labeyrie  Technique!  Anisoplanatism 


Our  starting  point  for  the  evaluation  of  the  effect  of  anisoplanatism 

on  the  Labeyrie  technique  is  Eq.  (60),  which  represents  the  spatial  frequency 
power  spectral  density  for  a pair  of  point  sources.  As  we  recall,  in 
Eq.  (60)  we  retained,  but  in  unreduced  form,  the  distinction  between 
effects  of  the  two  point  source  directions,  0 and  8 / , on  the  random 

J J 

phase  shifts,  0(0^;  x ) and  0 ( 0^'  ; x)  . At  this  point  we  return  to 
Eq.  (60)  but  with  the  objective  of  developing  statistical  results  showing 
the  dependence  of  the  power  spectral  density  on  the  angular  separation, 

-4  -4 

A _ 0 / 

'j  J 

We  shall  aim  our  algebraic  manipulation  of  the  terms  in  Eq.  (60) 
at  the  objective  of  obtaining  results  which  can  make  use  of  a statistical 
function  which  we  shall  call  the  hyper  wave- structure  function,  and 
which  we  define  as 

5>U,p)  = <[0(0+5;r+J)-0(0+5;r)] 

X [ 0 (0  ; ? + p ) - 0 (0  ; r ) ] > . (66) 

This  quantity  has  previously  been  studied  13  and  shown  to  be  expressible 
in  terms  of  an  integral  over  the  propagation  path.  It  can  be  written  as 

® ( «>  » P ) = 8*  16  (~zt~  ) S dv  C« ? f d*  C 1 " exP  ( 1 * * P ) 3 

' ' PATH 


X h"u/3  [ exp  ( i h • 5 V ) + exp  ( - in-  3 v ) ] , (67) 
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where  h is  a two-dimensional  (spatial  frequency)  vector,  and  v is  a 
variable  of  integration  measuring  position  along  the  propagation  path, 
with  v = 0 at  the  plane  of  the  optics  aperture,  where  we  "measure"  the 
wavefront  distortion  statistics.  C2  is  the  so  called  refractive -index 
structure  constant  and  is  a measure  of  the  local  optical  strength  of  the 
atmospheric  turbulence. 


-i<[0(0j;r)-0  ( eTj'  ; r')  ]s>  - I < [ 0 (0j  ; r ) - 0 (0/  ; r)  ]2) 

+ i < l 0 (0  j ; r ) - 0 ( 9j'  ; r + Xf  ) > + £ < [ 0 ( 0 } ; r ) - 0 ( 0 ;r  + Xf ) ]s  > 

+ i < L 0 (0,  ;r  + Xf  ) - 0 (?/  ; r')]2>  +|  <[  0 (0j  ,7  + Xf  ) - 0(?j/  ; r)]2> 

= < 0 ( 0J  ; r ) 0 ( 0 ,#  ; r* ) > + < 0 ( 0 J ; r ) 0 ( 0J  ; r ) > 

- I < 0 (0j  ; r ) 0 (?/  ; r'+  Xf ) > - \ < 0 (0J  ; r ' ) 0 ( 0,'  ; r + Xf ) > 

- i ; ? + X?)  0 (?,'  ; r')  ) - i < 0 (0J  ; r*  + X?)  0 («?j'  ; r ) > . (68) 


Our  problem  at  this  point  reduces  to  that  of  casting  the  six  terms 
still  explicitly  in  ensemble  average  form  in  Eq.  (60)  into  forms  that  can 
be  identified  with  Eq.  (66)  . The  procedure  is  tedius  but  because  of  its 
ultimate  importance,  is  worth  presenting  in  detail  here.  Starting  with 
the  expression  of  interest  in  Eq.  (60),  we  proceed  as  follows, 


In  obtaining  the  final  form  of  Eq.  (68)  we  have  made  use  of  the  fact  that 
the  wavefront  distortion  statistics  are  stationary  so  that  the  mean  square 
value  of  the  phase  shift  , ([0(0  ;x)]2),  is  the  same  no  matter  what 

the  values  of  * and  x . As  a consequence  all  of  the  mean  square  phase 
terms , when  summed  exactly  vanish.  Proceeding  in  the  same  way  we 
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can  show  that 


- < [ M9,  ; r ) - 0 (»;  ;•?')  ] [ 0 (e/;  r)  - 0 (e/;?')  3 > 

+ i < C 0 (03  J r)  - 0 (0j;  ? + xf  ) ] [ 0 (0^  ; r)  - 0 (0/  ; r + Xf  ) ] > 

+ i < [0(0J;r')-0(03  ; r + \i  ) ] [ 0 (8/ i?)  - 0 (?/ ; r + H?)  ] > 

= < 0 (93  ; r ) 0 (0/  ; r 7 ) > + ( 0 (0^  ; "?  ) 0 (7/  ; r ) > 

- i < 0 (0j  ; t ) 0 (0/;  T*7  + Xf  ) > - £ < 0 ( 03 ; 7 + Xf  ) 0 (03/ ; r ) > 

- |<  0 (0j  ; r' ) 0 (0  / ; r + Xf  ) ) - | < 0 (0  3 ; r + Xf  ) 0 (0/  ; ?' ) > . (69) 

In  developing  the  final  form  of  Eq.  (69)  we  have  again  made  use  of  the 
statistical  stationarity  of  the  random  wavefront  distortion.  This  allows 
us  to  argue  that  the  value  of  < 0 (0j  ; x ) 0 (0^  ; x ) ) is  independent 
of  3 . This  in  turn,  results  in  the  summing  to  zero  of  the  several  terms 
of  this  form  that  would  otherwise  appear  in  the  final  form  of  Eq.  (69)  . 

It  now  follows  directly  from  comparison  of  the  final  forms  of 
Eq.  's  (68)  and  (69)  that  the  two  starting  forms  are  equal.  The  starting 
form  of  Eq.  (68)  represents  a part  of  the  exponent  in  Eq.  (60),  while 
the  starting  form  of  Eq.  (69)  can  be  seen  , from  consideration  of  Eq.  (66), 
to  correspond  to  the  sum  of  several  hyper  wave- structure  functions. 

Thus  we  can  rewrite  Eq.  (60)  as 


</ (?)  = i | 31  | * JJdrdr7  W(  r ) W(  r+X?)  W(?7)W(-?'+  X?) 


X | C(A3S)a+  (A/*)3]  exp  (-{^r  - 7)  + *(X?)  -lUbG-r  -X?)t M?-?'  +X?)]]) 
+ 2(AJa)(AJ's)  cos  [2tt r?  • (0,-3/)]  exp  ( - (35(3, -<?/ , r-?')+J(\T) 


- * C3X?j-0/,  r-r' -X?)+S  (03  -0,/  , r'-7  - X?)]})| 


(70) 


Examination  of  Eq.  (67)  will  show  that  the  hyper  wave -structure  function 
possesses  sufficient  symmetry  that 

® (?,?)  = © U.  - p)  . (71) 

This  allows  us  to  rewrite  Eq.  (70)  in  the  form 

2(?)  = i 1 21 1 * J/dr  dr7  W(?)  W(?  + Xf)  W(r')  W(?  + X?) 

X | [(A,* ) + (A/a  )]2  exp(  - [Mr  - ?)+£  (X?)  [>(  r - r7-  X?)  +^(r  - r'+  X?)]  }) 

+ 2(A/)  (A3/a)  cos[2tt  ?•  (9  j - 0/  )] 
x exp^-flHej -03/ ,7-r7  ) +^(x?)  - iL®(0j  -0/,  X?) 


+ ®(0J  - 0,/  , r - r + Xf)]})| 


To  isolate  the  field-angle  dependence  we  rewrite  this  as 


(72) 


1 
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J (T)  = i | « | 4 J/dr  a?  W(r)  W(r  + \T)  W(r')  W(r'  + \T) 

X ||[(Aj3 )a  + (A/2)2]  exp(-U(r- "r')+i(U).^|r-?'.)i?)+^-r'+l?)]]) 


+ 2(Aj3)(Aj/3)  cos[2tt?  • (?3  -8j/  ) ] 

X expl-fXr-r^+^xfj-H^r-r'-Xfj+^r-r'  + Xf  )]3+Q(^  " */.  f-r',  Af  ))^»( 73 j 

where 

Q(0r0j/  , Xf ) = {Lfr(r-1')  - $ (8,  -0/,  ?-?' )] 

- -®(03-?3/,?-r'-X?)] 

- } Ulr-r'  + X?)  - S(03-  0j/  , ?-r'  + xl}]}.  (74) 

It  is  obvious  that  when  Q equals  zero  J 05  is  identical  to  what  we 
computed  in  the  previous  section,  when  there  was  no  anisoplanatism 
effect.  It  can  be  seen  from  a comparison  of  the  definitions  of  the  wave- 
structure  function,  Jfr,  and  the  hyper  wave- structure  function,  5),  as 

■4  ^ 

provided  by  Eq.  's  (19)  and  (66)  respectively,  that  when  6 } - 0 j*  is  equal 
to  zero  Q will  vanish  and  there  will  be  no  isoplanatism  problem.  The 
question  we  address  ourselves  to  is  how  large  - 9 }/  can  be  before  we 
encounter  anisoplanatism  problems,  i.  e.  , before  Q becomes  comparable 
to  or  larger  than  minus  one. 

We  recall  from  our  analysis  in  the  previous  two  sections  that 
when  we  are  interested  in  large  values  of  X|f  | , i.e.  , in  the  high  spatial 
frequencies  the  value  of  the  double  integral  defining  the  power  spectral 
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density,  J , is  determined  by  the  region  of  (r  , r* ) space  where 
| r " 1*  | i®  much  less  than  X|?  | . Accordingly,  we  restrict  our 
investigation  of  the  relationship  between  the  value  of  ? i - "e  / and 
anisoplanatism  to  the  case  where  X|f  | is  much  larger  than  r0  , 
i.e. , to  high  spatial  frequencies,  and  where  | r - "rf  | is  small,  in 
particular,  much  smaller  than  X \t  | . Making  use  of  Eq.  (67),  hyper 
wave- structure  function,  S , and  noting  that  when  3 s 0 it  also 
provides  a definition  for  the  ordinary  wave -structure  function,  jfr  , 
we  can  write  Q as 


Q («? » "p  » X? ) = ^ k®  J* dv  CNS  J*d * k “ii/3 

4TT  PATH 


{^exp[iH  • ("p  + X?)]  - exp[in  • p]  + £exp[in  • - X?)]] 


X {-exp(i>t  • tfv)  + 2 - exp(-in  • ^v)] 


= " k*  J* dv  Cn2  Jd«  h"11/3  exp(in  • p ) 

path 

X [exp  (iK  • \?/2)  - exp  (-in  • Xf  /Z)f 
X [ exp  (in  • tfv/2)  - exp  (-in  • ?v/2)]S 
= kS  J*  dvCNa  Jd^'n/3  exp(i*.7) 

T PATH 

X sin3  \?  • h)  [ exp  <i >t  • 3v/2)  - exp  (-  in  • Jv/2)]3  . (75) 


At  this  point  we  note  that  X(f  | is  very  large  compared  to  r0  and  thus 
large  with  respect  to  all  of  the  values  of  p and  of  tfv  of  interest.  As  a 
consequence  we  may  argue  that  the  sin3  ( $ X?  • k ) term  goes  through 
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many  oscillations  in  any  locality  of  the  x-  space  integration,  and  there- 
fore may  be  considered  to  have  an  approximate  value  of  one-half.  Thus 
we  may  rewrite  Eq.  (75)  as 

Q Q,  "p,  Xf ) = k3  S dv  CH  I d*  H"u/3  exp  (in  • p ) 

PATH 

X [exp  (in  • £v/2)  - exp  (-in  • £ v/2)]2  . (76) 

We  can  rewrite  this  as 

Q (£»  "p,  X?)  = - "g—  k2  j*  dv  CN2  J d"x  x”11/3  cos  ( n • "p) 

71  PATH 

X sin2  (x  • tfv/2)  , (77) 

where  we  have  dropped  the  imaginary  part  of  exp  (in  • 'p)  since  it  is 
an  odd-function  of  x while  all  the  rest  of  the  integral  is  an  even-function 
of  x , so  that  over  the  infinite  range  of  the  x -integration  its  contribution 
would  vanish. 

We  now  further  argue  that  the  constraint  that  we  are  only  interested 
in  small  values  of  |"p  | , since  only  small  values  of  |r  - r*  | contribute 
significantly  to  the  value  of  J , allows  us  to  approximate  the  isoplanatism 
constraint  by  requiring  that  Q ($)  , where 

G (»?)  = k2  J dv  CN2  ^dx  x “u/3  sin2(x  • #v/2)  (78) 

17  ~ PATH 

be  less  than  unity.  Here  we  have  replaced  cos  (x  * "p)  with  unity  based 
on  the  fact  that  with  |"p  | small  or  comparable  to  rQ  , x • "p  will  not  be 
large  for  any  value  of  x which  contributes  significantly  to  the  wavefront 
distortion. 


Letting  v denote  the  angle  between  >?  and  £ we  can  rewrite 
Eq.  (78)  as 


^ W = n2  k2  Jdv  CNS  J*  d*  h"^3  J dv  sin2  [£  v cos  (v)] 


2TT 


PATH 


0 

an 


8 16  2 00 

2na  k Jdv  CN2  Jd*  J dv  { l-cos[ntf  v cos  (v)]} 

PATH  0 0 


k2  ; dv  CN2  JdK  H-a'3  [1  - J0(k*v)] 

PATH  0 


TT 


(79) 


Making  a change  of  variables,  with  x = x«?v,  we  can  rewrite  Eq.  (79) 


as 


Q (tf)  = d6/3  k2  J dvyS/3  CN8  Jdxx"®/3  [1  - J0  (x)]|  . (80) 

It  is  easy  to  show  1 4 that 


Jdx  x-*»  [1  - J0(x)]  = - Z-a'3 


where 


= 4 2 


dv  v6^  C 


We  call  t?0  the  isoplanatism  angle. 

In  Fig.  6,  we  show  a sample  vertical  distribution  of  the  optical 
strength  of  turbulence,  CNa  , as  reported  by  Greenwood.  16  Using  this 
distribution,  we  have  calculated  the  isoplanatism  angle  for  a wavelength 
\ = 5.  5 X 1CH’  m . We  find  that  in  this  case  t?0  = 8.  6 X 10*«  rad  . Our 
theory  predicts  that  the  Labeyrie  technique  should  produce  good  results 
when  working  with  a pair  of  point  sources  with  angular  separation  le3S 
than  this,  but  that  if  the  point  sources  are  separated  by  an  angle  greater 
than  this,  the  high  spatial  frequency  details  will  be  lost,  or  at  least  sig- 
nificantly attenuated  from  what  we  expect,  based  on  calibration  of  the 
technique  with  a single  point  source.  Referring  back  to  Eq.  (73)  and 
noting  that  Q is  the  negative  of  Q , we  would  expect  from  considera- 
tion of  Eq.  's  (82)  and  (84)  that  the  high  frequency  details  will  be  attenu- 
ated approximately  as  exp  [ - (j>  /#0)&/i3  ] relative  to  the  ability  of  the 
Labeyrie  technique  to  accommodate  a single  point  source.  From  con- 
sideration of  Eq.  (73),  we  can  see  that  the  net  effect  of  a lack  of  isoplana- 
tism will  be  to  make  the  source  appear  more  like  a single  point  source. 

A pair  of  equal  intensity  point  sources  would  appear  like  a pair  of  un- 
equal intensity  point  sources,  the  ratio  of  intensities  appearing  to  be 
£ exp  [-G»/0o)b/S]  . 

With  these  results  in  hand,  we  are  now  ready  to  turn  our  attention 
to  the  Knox-Thompson  concept.  In  the  next  section,  we  shall  describe 
this  concept  and  then  present  an  analysis  of  its  operation,  assuming  that 
we  are  dealing  with  a pair  of  point  sources  close  enough  together  that 
there  is  no  anisoplanatism  effect.  Then  in  Section  6,  we  shall  extend  this 
analysis  to  treat  the  case  where  the  angular  separation  is  large  enough 
for  anisoplanatism  effects  to  occur. 
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Knox-Thompson  Algorithm:  Isoplanatism  Assumed 
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In  the  absence  of  anisoplanatism  effects  and  with  a "calibration 
run"  against  a known  point  source  the  Labeyrie  technique  will  produce 
a useful  and  valid  spatial  frequency  power  spectrum  of  an  object  of 
interest.  While  this  gives  us  much  of  the  information  we  need  to  know  about 
the  object,  it  does  not  permit  formation  of  the  image  of  the  object.  We 
have  information  on  the  amplitude  of  each  spatial  frequency  of  the  image, 
but  we  have  no  knowledge  of  the  spatial  phase  shift  to  be  associated  with 
each  spatial  frequency  — so  we  cannot  reconstruct  the  image. 

The  basic  problem  in  determining  the  phase  shift  stems  from 
the  fact  that,  except  for  the  very  lowest  spatial  frequencies,  (namely 
|?|  < r0/'X  ) , the  turbulence  induced  part  of  the  phase  shift  introduces 
an  rms  spatial  phase  shift,  a , that  is  many  radians.  As  a consequence 
the  average  of  the  random  spatial  frequency  component,  S(f ) , will  be 
so  strongly  attenuated  that  we  cannot  use  the  ensemble  average, 

(S  (f)  ) , to  determine  the  part  of  the  phase  shift  due  to  the  image.  The 
attenuation  is  exp  (-a2)  , which  is  too  small  to  allow  any  practical 
use  to  be  obtained  from  the  ensemble  average  ( S (?)  ) . 

In  mathematical  terms  we  would  write  that 

S (?)  = A(At  exp  [i  (cp,  + <PT)]  , (85) 

where  A,  is  the  image  associated  amplitude,  (which  we  can  determine 
using  the  Labeyrie  technique),  and  cp ( is  the  image  associated  phase  shift 
which  we  wish  to  determine.  AT  is  the  random  amplitude  factor 
introduced  by  turbulence  effects  and  cp  T is  the  turbulence  induced  spatial 
phase  shift,  whose  variance  is  of  , If  we  determine  the  (ensemble) 
average  of  the  random  amplitude,  S (f  ) , we  get 
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<S(f»  = C < AT  > exP  (-CTfa)  } A,  eacp  (i  <p,  ) . 


(86) 


L 


The  spatial  phase  shift  information  is  there  and  we  could  easily  extract 
it  as  the  argument  of  the  complex  value  of  ( S (?)  ) , but  because  of  the 
exp  (-c  2 ) factor  we  are  dealing  with  such  an  unfavorable  signal  strength 
situation  that  we  cannot  obtain  a useful  estimate  of  the  spatial  phase 
shift.  Our  answer  will  be  determined  by  the  noise  rather  than  by  the 
signal. 

An  alternate  approach  for  the  calculation  of  the  image  associated 
spatial  phase  shift,  <p,  , has  been  suggested  by  McGlamery16.  This 
involves  extracting  the  phase,  (cp  , + <PT  ) , from  each  sample  of  the 
random  spatial  frequency  complex  amplitude,  S (? ) , and  then  (ensemble) 
averaging  this  phase  shift  to  obtain  < ( <p  , + <pT  ) ) which,  since  <<PT  ) = 0 , 
should  reduce  to  ((<P,  +CPT))=<P,  . Unfortunately,  due  to  a combin- 

ation of  noise  effects,  [irrevocable  incorporated  into  the  results  by  the 
nonlinearity  of  the  process  of  calculating  the  argument  of  S (?)  ] , and 
due  to  a basic  2n- ambiguity , there  is  an  uncertainty  of  the  order  of 
2tt  N in  our  knowledge  of  (9,  + <PT  ) , and  as  a consequence  we  do  not  get 
a meaningful  result  when  we  attempt  to  detei-mine  ( (<P,  + «PT  ) ) . 

The  calculation  of  the  spatial  frequency  power  spectral  density 

•f  (?)  = < S*  (?)  S (?)  > , (87) 

does  not  encounter  any  difficulty  due  to  the  turbulence  induced  random 
phase  shift,  cpT  , since  there  is  a cancellation  of  this  term  in  taking  the 
product  S*(f)S(f).  The  Knox-Thompson  concept  seeks  to  get  around 
all  of  the  above  difficulties  by  measuring  the  bispectrum 
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. ?')  = <S*(?)S(?)  > , 


(88) 


I 


where  ? and  1'  are  distinct  but  not  very  widely  different  spatial 
frequencies.  As  will  be  seen  in  this  case,  the  two  turbulence  induced 
phase  shifts,  <pT  and  cpT'  , are  closely  correlated.  From  this  it 
follows  that  <PT  - <p/  > which  is  the  turbulence  induced  random  argument 
of  the  complex  product  S*  (?'  ) S(?  ) , has  a small  mean  square  value 
aff • , and  so  the  attenution  factor,  exp  (-o3  t*  ) , encountered  in  forming 
the  (ensemble)  average  is  of  the  order  of  unity.  The  condition  required 
for  this  to  be  the  case,  i.e.  , for  <jf  f # to  be  very  small,  will  be  seen  to 
be  |?  - T'\  < r0/\  . 


We  shall  start  our  analysis  of  the  Knox-Thompson  concept  by  con- 
sidering a pair  of  point  sources  with  amplitudes  and  A ^ , located 

at  angular  positions  9 and  0 . Making  use  of  the  intermediate  re- 

suits  of  Section  3,  in  particular  of  Eq.  (52),  we  see  that  we  can  write  the 
bispectrum,  as  defined  by  Eq.  (7),  as 

;f,g)  = 1 |«|*  JJ  dr  d?'  W(r)  W(r+\?)  W(?')  W(r'+x?') 

X ( •|(AJ8)*  exp  [ -2n  i (f -f  n)  exp  {i[0(0,j;"r)  - 0 (<f3  ;"?+ >T) 

+ 0^;t'+\7')  -0^;r>)]} 

+ A3a  A*  exp  [ - 2tt  i (f  • ^ -7*^)3  exp  {<t>($}\r) 

-0(?3;r +\T)  +0(fJ,;r'+xf'')  -0(^,;r')]} 

+ A33  Aj3  exp  [-2tt  i ~ '9})1  exP  (i[0(^:'?) 

- 0(?3,;r+xf')  + 0(f3  ;r'+\?')  - 0(?3  ;r')]) 

+ (A  33)*  exp  [-2tt  i (T-T')  exp  (i  [0  (fTj,;"r)  - 0(ej,;‘r+\T') 

-0  (^r;r'+x?>>)  - 0(e3,;r')]}|.> 
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(89) 


Here  we  have  written  a product  of  integrals  as  a double  integral,  carried 
out  the  multiplication,  and  grouped  terms  as  appropriate  in  the  integrand. 

We  now  introduce  the  assumption  that  the  statistics  of  0 CTj ; r ) , 
with  respect  to  0($j,;'r>)  are  identical  to  the  statistics  of  r)  with 

respect  to  0(t;  r ')  — in  essence,  assuming  that  and  "9  are  in  the 

3 3 3 

same  isoplanatic  patch.  Thus  restricting  our  attention  in  the  balance  of 
this  section  to  the  case  where  isoplanatism  is  assumed,  and  making  use 
of  Eq.  's  (16),  (17),  and  (19),  we  can  rewrite  Eq.  (89)  in  the  form 

7')  = % |a|*  ff  d?  dr  ' W (?)  W(r+\7)  W(?')  W(?'+\T') 
x exp  .|-£  [M\?)  + MX?')  +£(r-r') 

+#(?-? '+\<f-F')]  -.*(?-?'+\n  - ^(r-r'-xT')}  J. 

X ((A»)»  exp[-2Tri  (?-?')  •?,] 

+ {Aa)(Af»)  exp  [-2rri  (T*^  -T'*^,)] 

+ (A  *)(A}»)  exp  [-2n  i - T-fy] 

+ exp  [ -2tt  i (F-?')  •*,,]}  . (90) 

We  can  rewrite  this  as 

= A 1 9||4  fj  dr  dr'  W(r)  W(r+\T)  W(?')  W(?'+\?') 

X exp  {£{x?)  + + £(?-"?')  + '+\(T-?')] 

- - ;>(?-?'  - xT')}J 

X {[(A*)  exP  (-2ni  + (Aj*)  exp  (-2tt  i?  •tfj,] 

X [(A»)  exp  + (Aya)  exp  (2tt  iT*  *1^)] ) . (91) 
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It  is  particularly  significant  to  note  here  that  the  two  square  bracket  terms 
in  the  final  curly  brackets  of  Eq.  (91)  represent  the  fourier  transform  of 
the  two  point  source  pattern  at  spatial  frequency  f (in  the  first  of  the  square 
brackets),  times  the  complex  conjugate  of  the  fourier  transform  of  the  pat- 
tern at  spatial  frequency  7'  (in  the  second  square  bracket).  As  noted  pre- 
viously, the  phase  shift  associated  with  this  product  is  the  difference  of  the 
phase  of  the  two  spatial  frequency  components,  at  7 and  7 * , of  the  object 
pattern  — and  if  we  can  determine  this  difference  of  phase  shift  for  all 
"adjacent"  pairs  of  spatial  frequencies  in  a matrix  covering  spatial  frequency 
space,  we  can  determine  the  absolute  phase  shift  for  each  spatial  frequency 
component.  (The  process  is  essentially  equivalent  to  that  described  by  us 
elsewhere17 for  determination  of  wavefront  distortion  from  an  array  of  wave- 
front  tilt  measurements.  ) 

The  key  question  is  whether  or  not  we  can  determine  the  difference 
of  phase  shifts  from  a measurement  of  the  bispectrum,  »?(f  ,7')  . This, 
in  turn,  is  determined  by  the  value  of  the  "r,  7 '-double  integral  in  Eq.  (91). 
Comparing  the  integrand  here  with  the  corresponding  integrand  in  Eq.  (63), 
for  the  Labeyrie  speckle  interferometry  technique,  we  see  that  the  two  will 
be  equal  when  7 equals  7 . If  T'  is  nearly  equal  to  7 , then  the  inte- 

grals should  have  nearly  the  same  value  — but  exactly  what  difference  be- 
tween r and  7*  is  allowable  under  the  expression  "nearly  equal.  " We 
recall  that  in  the  series  of  steps  leading  up  to  Eq.  (40),  and  thus  also  to 
Eq.  (63),  it  was  shown  that  when  7 and  7'  are  exactly  equal,  the  integral 
had  a value  of  0.435  (rQ/D)»  t (f)  . In  the  following,  we  shall  repeat  those 
steps  with  7*  nearly  but  not  exactly  equal  to  7 , and  show  that  if  \|f  -7'  | 

« rQ  , then  the  value  of  the  double  integral  will  be  essentially  as  large  as 
when  7*  is  exactly  equal  to  7 . 

It  is  convenient  to  start  by  defining  the  quantity  g by  the  expression 
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<f  = 'i  UOI)  + M\T')  + + rfr-t'+xf-T')) 

Making  use  of  Eq.  (22),  we  can  write 

* * 3-  44  r0-BA»[xf)B>*+  (xf')B*  + (■?-?' +x(T-F')|«^ 

-!?-■?'+  xf|**  - (?-■?' -xr'p*]  , (93) 

which  expression  we  may  usefully  compare  with  Eq.  (23).  In  our  previous 
treatment  of  the  Labeyrie  technique,  developed  from  Eq.  (23),  with  Eq.  (18) 
playing  the  same  role  as  Eq.  (21)  in  defining  the  range  of  values  of  ? and 
■?'  to  be  considered  and  the  weighting  to  be  assigned  to  each  (?,"?') -value, 
we  found  that  when  Xf  was  much  greater  than  rQ  , then  because  there 
were  some  very  "convenient"  cancellation  of  terms  in  the  approximation  for 
Eq.  (23)  when  |r-‘?-'|  was  much  less  than  Xf  . it  followed  that  the  value 
of  the  (r,’?')- integration  was  determined  by  the  range  of  integration  in  which 
|r-"?  | was  less  than  or  about  equal  to  rQ  . In  the  same  way,  we  may  ex- 
pect here  that  when  there  is  negligible  difference  between  the  values  of 
T and  T'  , and  is  much  larger  than  rQ  , then  the  value  of  the  (?,"?')- 
integration  in  Eq.  (91)  will  be  determined  by  the  contribution  from  the  region 
in  which  I?--?')  is  less  than  or  about  equal  to  rQ  . This  means  that  as  a 
practical  matter  in  evaluating  the  bispectrum,  for  \i  much  greater 

than  rQ  , and  T and  T'  nearly  equal,  we  can  approximate  $ from  Eq.  (93) 
subject  to  the  assumption  that  is  less  than  or  about  equal  to  r^  . 

Making  use  of  Eq.  (25)  and  equivalent  approximations,  we  can  now 
rewrite  Eq.  (93)  as 

<#  = 3.  44  r0-e*  ["  + |M'+x<T-r')|«* 

+ | (xf H*(?--?') • \(T -T')  + | (XfH*(?--?') • \T  . (94) 
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It  is  immediately  apparent  from  consideration  of  Eq.  (94)  substituted  into 
Eq.  (91)  that  the  key  question  in  determining  the  value  of  the  integral  is 
whether  X ( f is  greater  than  rQ  , or  less.  In  either  case,  the  result 
is  not  significantly  dependent  on  anything  except  the  |‘r-r'|**i  and  the 
)‘?-7#+  \(f-f')|*/!s  terms  in  Eq.  (94).  We  can  rewrite  Eq.  (91)  as 

J&.7')  |»j|*  Jf  dr  <&'  W(?)  W(?+xf)  W(?')  W(f'+\7 ') 

X exp  {-3.44  r0-e/3[  |r-7'|6/3+  |7-?'+  X^-^)|e/3]  j 
X {[(A  ■)  exp  (-2rr  iF*^)  + (A}»)  exp  (-2tt  if  • (J^)] 

X [(A}8)  exp  (2tt  if' -1^)  + (Aj8)  exp  (2tt  iT' . (95) 

For  the  required  condition  for  successful  operation  of  the  Knox-Thompson 
technique,  i.  e.  , for  X(f-f')  much  less  than  rQ  , and  recalling  that  the 
value  of  the  integral  is  determined,  for  the  most  part,  by  the  region  in 
which  |?-7'|  is  of  the  order  of  rQ  , we  see  that  Eq.  (95)  reduces  to 

j(f,7')**\  |Hl4  JJ*  dr  d?'  W(?)  WCr+xf)  W(?')  W(r'  + \7') 

X exp  [-6.88  (|r-'r'| /r0)6^] 

X [[(Aj8)  exp  (-2rr  if  *^)  + (Aj3)  exp  (-2tt  if -f^)] 

X [ (Aj3)  exp  (Zvil*  -fj)  + (Aj8)  exp  (2nif  *7^)]}  . (96) 

This  integral  was  previously  encountered  in  Eq.  (30)  and  evaluated  in 
Eq.  (39).  Making  use  of  Eq.  (21),  we  see  that  with  this  integral  evaluation, 
we  can  write  the  normalized  bispectrum  as 

= t°-435  ,ro/D)“  V<r» 

X {[(A3)  exp  (-2ni?*fJ)  + (A^3)  exp  (-2tt  if  • fJ#)] 

X [(Aj3)  exp  (2tt  if'  -fj)  + (A^8)  exp  (2n  if'  ’l^)] 

X C(Aj8)  + (Aj8)]-3  } . (97) 
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The  leading  curly  bracket  term  is  the  quantity  of  interest.  It  represents 
the  ability  of  the  process  to  "avoid"  washing  out  the  information  carrying 
part  of  the  signal  because  of  averaging  over  a large  range  of  phase  fluc- 
tuations. The  larger  the  curly  bracket  term  is,  the  better  the  system  may 
be  considered  to  perform.  We  note  that  with  X|F-P|  smaller  than  r0  , 
the  system's  ability  to  avoid  signal  washout  is  as  good  as  for  the  Labeyrie 
speckle  interferometry  technique. 

If  we  violate  the  requirement  that  X|F-F'|  is  to  be  less  than  r0  , 
then  since  the  value  of  the  double  integral  in  Eq.  (95)  comes  principally 
from  the  (?,?*)  region  of  integration  where  j is  less  than  r0  , 

we  can  approximate  Eq.  (95)  by  the  expression 

^r.F-)  = exp  [-3.44  (x|7-r'|/re)6/3]  j |V|4 

X J*J  dr  dr'  exp  [-3.  44  ( |"r-"r  ' | /r 0)P^ 

X W(?)  W(?+xF)  W(?')  W(r'+XF') 
x [[(A,3)  exp  (-2rriF*^)  + (A^3)  exp  (-2rt  iF*^,)] 

X [ (Aj3)  exp  (2tt  \T*  ‘fy  + (Aji)  exp  (2tt  iF'  •FJ,)]  ) . (98) 

Except  for  a slight  change  of  coefficient  in  the  exponent  [which  we  can 
"get  rid  of"  by  replacing  3.44  r0-«*  by  6.88  (1.  516  r0)-«*  ] , the  integral 
is  the  same  one  which  we  just  evaluated.  Carrying  out  the  integration  and 
normalizing  as  before,  in  this  case  we  get 

X {[(Aj3)  exp  (-2tt  iF*^)  + (A^ •)  exp  (-2ni F*-^)] 
x [ (A^3)  exp  (2tt  if'  *7^  + (A»)  exp  (2tt  if'  *1^,)] 

X [(A,3)  + (A^3)]-3  ) (99, 


l 


- 42  - 


It  is  immediately  apparent  from  inspection  of  this  result  that  if  we  allow 
the  frequency  difference  to  be  large,  i.  e.  , if  \|fr-?'|  is  larger  than  rQ  , 
the  Knox-Thompson  technique  will  produce  poor  results,  i,  e.  , the  inform- 
ation carrying  part  of  the  signal  will  be  washed  out. 

From  a comparison  of  Eq.  's  (97)  and  (99),  we  can  form  an  estimate 
of  how  much  smaller  than  rQ  we  should  require  Xl?-?')  to  be.  We  form 
our  estimate  based  on  the  requirement  that  the  knee  of  the  curve  of  (7-7* J 
dependence  should  correspond  to  the  value  where  the  two  (asymptotic)  depen- 
dencies intersect.  This  occurs  where 

exp  [-3.44  (X | f-F'l/r^/a  = 0.435  , (100) 

where  the  subscript  K denotes  the  "knee"  value.  Solving  this  equation, 
we  get 

X|7-f'|K  = 0.427  r0  . (101) 

Based  on  this  result,  we  would  suggest  that  for  proper  implementation  of 
the  Knox-Thompson  algorithm,  the  data  processing  should  be  done  with  an 
array  of  spatial  frequencies  whose  frequency  spacing  is  no  greater  than 
0.4  i..  ’ , and  preferably  no  greater  than  0.2  rQ/X  . 

With  this  result  in  hand,  we  are  now  ready  to  proceed  to  the  analysis 
of  the  Knox-Thompson  algorithm  in  the  case  where  the  angular  separation, 

-Ij,  , of  the  two  point  sources  is  large  enough  that  there  will  be  aniso- 
planatism  effects.  We  take  this  up  next,  in  Section  6. 


6.  Knox- Thompson  Algorithm:  Anisoplanatigm 

We  start  our  analysis  of  the  effect  of  anisoplanatism  on  the  Knox- 
Thompson  algorithm  with  Eq.  (89).  Here,  unlike  our  procedure  in  the 
last  section,  we  shall  not  introduce  the  assumption  that  the  wavefront  dis- 
tortion, 0{?fr)  , is  independent  of  the  exact  value  of  F . Making  use  of 
Eq. 's  (16),  (17),  and  (19),  we  can  recast  this  as 

jtf.T')  = $ |a|*  .fj  d?  dr'  W(r)  W(?+X?)  W(?')  W(?')  W(?'4.XF') 

X |(A»)»  exp  [-2*1  (T-F')-^]  exp  U(\F)  - fitf-V'-tf') 

+ ^(?-r')+ ^r--?'+X(r-r')]  -^(?-?'+x?)+^(Xr')}) 

+ A»AJ»  exp  C-2ni(r-0J-r'-fJJ]  exp  (-£  U(X?)-<[0(e3;?) 

- 0(?r;?'+xr')]s»  > 

+ <[0(F3  ,-r)  -0^;7')]3  )4<[0(^  ;r+XF)-0  '4X?')]»  ) 

-<[0(9j  ;r+X?-0(?j,;r')]»  > -I-  £()J')}^ 

+ AJ*Aj*  exP[-2"itf'-Vr'-?j)]  exp  (-^U(xr)-<[0(yj;?) 

- 0(?J;'?'+X?']8) 

+ a<t>(frir)-<p(S}  ;■?')]•  >+<[0^;?+xf)-0(S;;?'+xr')]»> 

- <[0(^,/.r+xF)-0(^  ;"?')]*  ) + ^(xf')}^ 

+ (a3=)»  exp  [-2ni(r-r').*r]  exp  (-iU(xr)-^{7-F'-Xr') 

+^(?--?')+^[-?-7'+x(r-?')]-^(?-?'+xr)+^(X?'))N}  . 

We  now  make  use  of  the  symmetrization  argument  used  to  let  us  go  from 
Eq.  (58)  to  Eq.  (60).  We  then  note  that  we  can  now  repeat  essentially  the 
same  manipulations  we  used  in  Eq.  (68)  and  Eq.  (69).  The  equivalence  of 


(102) 
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these  equations  allows  us  to  rewrite  our  results  in  terms  of  the  hyper 
wave -structure  function,  as  defined  by  Eq.  (66),  so  that  we  obtain  in 
place  of  Eq.  (102)  the  result  that 

\ |M|*  J*J*  dr  d?'  W(?)  W(?+xD  W(?')  W(?'+X?') 

X {[(A»)»  exp  [-2ni^-?')*l;]+(AJ*)*eXp[-2TTia,.?')*l;j] 

x e*P  (-\  {^(Xn+^(X?')+^(r-r ')+^[r-?'+\(T-P)] 

-^(r - r '+  X?")  -^(r  - r ' - xf' ) } ^ 

+ (AJ»)(AJ  •) [exp  [-2ni(?-fi  -7'  • jj J]  + exp  [ -2n  i (F • -F'  .^ )] } 

x exP  (-i  (.*(x?)+MxF' -i;,;?-r')+5>[^  -■?J,;-?--?'+x(?-f')] 

-^-^.•^-r'+X?)  -®(^-e3,;?-r'-x?')}^j.  . (103) 

It  is  useful  to  rewrite  this  equation  as 

•rtf,**)  = i I Ml*  J*J*  d?  dr'  W(r)  W(?+xF)  W(?')  W(F'+xF') 


x exp  {^(Xf*)+^(xr')+^f?-?')+^[?-?'-|-xr-f#)] 

-^(r-"r'+xT)  -^(r-?'-X?#)}  ^ 

X ((A  ■)•  exp  [ -2n  i (T -7') • ^ -1JJ  + (As«)8  exp  [-2tt  i (F-F')  •*  ] 

+ (Aj8) (Aj  •)  {exp  [ -2n  i (?  • ^ -P  • J]+  exp  [ -2tt  i (T  • ^ -F  • ^ )]  } 
x exp  [in  (^  -^.r-r',  XF,  X?')]}^  , (104) 

where 

0(3;  -l;  XF,  \F')  = -srt; 

+ UCr-?'+x(F-r')]  - -^,;r-?'+X(F-F')]} 

- [/>(?-r'+>f)  - -pdj  -9r;r-r'  + xJ)] 

- tj>(r-?'-\7')  - -fr -\T#)]^  . (105) 


A comparison  of  Eq.  (104)  with  Eq.  (90)  makes  it  clear  that  the  condition 
on  |fa  | re<luired  in  order  for  there  to  be  no  significant  anisoplanatism 
effect  is  that  $ must  have  a value  close  to  zero  in  the  relevant  region  of 
the  (?,"?')  integral. 

To  establish  the  allowable  range  of  values  of  -‘j’  J for  isoplana- 
tism  to  apply,  we  shall  carry  out  an  evaluation  of  ofif.'p,  \T,  \T')  subject 
to  the  conditions  that  p is  less  than  or  of  the  order  of  rQ  , that  Xf  and 
Xf'  are  both  much  greater  than  rQ  , and  that  X|f-f'|  is  less  than  rQ 
— a 8 these  are  the  conditions  we  developed  in  the  last  section  as  defining 
the  range  in  f,7'  and  (r, ~r* -integration  space  within  which  we  get  signi- 
ficant contributions  to  the  bispectrum,  ,f')  . 

Making  use  of  Eq.  (67),  we  can  write 


©(*.  p.  Xf , XT')  = ^2.  » J dv  c*  J d£  H-u* 

P AT  H 

X (exp  [i»t • (p  + \?)3  + exp  [in  • Cp  - xf ')] 

- exp  (in-'p)  - exp  [in  • ["p + X(f-f')]}^ 

X {-exp  (i*  *#v)  + 2 - exp  (-ijt-^v)} 

With  a bit  of  trigonometric  manipulation,  this  can  be  rewritten  as 

o(*.  P.  xf.  xf')  = - k»  f dv  CNa  J <5*  h-u/3 

PATH 

X exp  {in  ’ x(f-f')]} 

X {sin*[iJ,H  • x(f-f')]  - sina[iH*x(T+r')]} 
X [exp  (in*?  v/2)  - exp  (-in’?  v/2)]8 

We  may  now  proceed  almost  exactly  as  we  did  following  Eq.  (75). 


(106) 


(107) 
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We  start  by  noting  thr.t  f or  virtually  all  of  the  jt-valuee  that  con- 
tribute significantly  to  the  wavefront  distortion  process,  £ \.(f4?')  is  so 


large  that  the  sirf  [£"k  • X(T+?')]  term  will  oscillate  rapidly  and  therefore 
will  have  an  average  value  of  one-half.  Thus  we  can  replace  that  sinc- 
function  squared  by  one -half.  Also,  we  note  that  since  all  of  the  rest  of 
the  integrand  is  manifestly  symmetric  in  h , we  can  replace  exp  {in* 

["p  + ^ X^+T")]  } by  the  corresponding  cosine -function.  Thus  we  can  write 
in  place  of  Eq.  (107) 

nfr,  p,  xf . )?')  = 1J^p-  ^ J dv  CNa  J <5?  h-m/* 

f AT  H 

x cos  sin^^K-xr-f')]} 

X sina(H*7v/2)  . (108) 

From  our  previous  work,  we  know  that  we  are  only  interested  in  values  of 
I"?-?']  , and  thus  of  "p  which  are  less  than  or  equal  to  rQ  . Similarly, 

we  know  that  we  must  restrict  our  attention  to  the  case  where  \|?-I>>|  is 
less  than  rQ  . On  the  other  hand,  we  know  that  only  values  of  h which 
are  smaller  than  rQ~*  contribute  significantly  to  the  wavefront  distortion. 

It  thus  follows  that  to  the  level  of  a first  approximation,  we  may  treat  the 
cosine-term  and  the  one-minus-two-times-sine-squared-term  in  the  inte- 
grand of  Eq.  (108)  as  though  they  had  values  about  equal  to  unity. 

Taking  note  of  these  approximations  [and  recalling  the  factor  of 
one -half  in  the  final  exponent  of  Eq.  (104)],  we  see  that  the  condition  on 
£ for  there  to  be  no  significant  anisoplanatism  effect  is  that  the  function 
O(0)  , where 

0(0)  = k»  f dvCN»  J dHK-u*  sin»(>?  v/2)  , (109) 

PATH 


a Vi  rail  have  a value  than  two.  Maltino  reference  to  Ea  . f 78 1 and  the 


discussion  there,  we  see  that  we  have  arrived  at  exactly-  the  same  condi- 
tion on  £ , or  |^-^|  * for  isoplanatism  to  apply  for  the  Knox-Thompson 

algorithm  as  we  previously  developed  for  isoplanatism  to  apply  to  the 
Labeyrie  technique.  We  require  that  Eq.  (83)  with  Eq.  (84)  be  satisfied 
for  the  Knox-Thompson  algorithm  to  yield  isoplanatic  performance. 
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Figure  1.  Speckle  Transfer  Function  for  D/r0  = 10  . 

The  rapidly  decreasing  curve  represents  the  low  frequency  approx- 
imation of  Eq.  (44).  The  everywhere  relatively  low  level  curve  represents 
the  high  frequency  approximation  of  Eq.  (40).  The  curve  running  the  full 
range  of  spatial  frequencies  represents  the  Monte  Carlo  results.  The 
composite  of  the  two  approximations  can  be  seen  to  constitute  a fair  approx 
imation  to  the  full  range  result  obtained  by  Monte  Carlo  methods. 
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Figure  2.  Speckle  Transfer  Function  for  D/r0  = 30  . 

The  rapidly  decreasing  curve  represents  the  low  frequency  approx- 
imation of  Eq.  (44).  The  everywhere  relatively  low  level  curve  represents 
the  high  frequency  approximation  of  Eq.  (40)  . The  curve  running  the  full 
range  of  spatial  frequencies  represents  the  Monte  Carlo  results.  The 
composite  of  the  two  approximations  can  be  seen  to  constitute  a fair 
approximation  to  the  full  range  result  obtained  by  Monte  Carlo  methods. 


SPECKLE  TRANSFER  FUNCTION 


Figure  3.  Speckle  Transfer  Function  for  D/r0  = 100  . 

The  rapidly  decreasing  curve  represents  the  low  frequency 
approximation  of  Eq.  (44).  The  everywhere  relatively  low  level  curve 
represents  the  high  frequency  approximation  of  Eq.  (40).  The  curve 
running  the  full  range  of  spatial  frequencies  represents  the  Monte  Carlo 
results.  The  composite  of  the  two  approximations  can  be  seen  to  constitute 
a fair  approximation  to  the  full  range  result  obtained  by  Monte  Carlo  methods. 
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Figure  4.  Speckle  Transfer  Function  for  D/r0  = 


The  rapidly  decreasing  curve  represents  the  low  frequency  ap 
proximation  of  Eq.  (44).  The  everywhere  relatively  low  level  curve 
represents  the  high  frequency  approximation  of  Eq.  (40).  The  curve 
running  the  full  range  of  spatial  frequencies  represents  the  Monte 
Carlo  results.  The  composite  of  the  two  approximations  can  be  seen 
to  constitute  a fair  approximation  to  the  full  range  result  obtained  by 
Monte  Carlo  methods. 
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Figures.  Speckle  Transfer  Function  for  D/r0  = 


■ The  rapidly  decreasing  curve  represents  the  low  frequency  ap 
proximation  of  Eq.  (44).  The  everywhere  relatively  low  level  curve 
represents  the  high  frequency  approximation  of  Eq.  (40).  The  curve 
running  the  full  range  of  spatial  frequencies  represents  the  Monte 
Carlo  results.  The  composite  of  the  two  approximations  can  be  seen 
to  constitute  a fair  approximation  to  the  fuli  range  result  obtained  by 
Monte  Carlo  methods. 
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Appendix  A 

Monte  Carlo  Integral  Evaluation 

In  order  to  develop  a set  of  values  for  the  frequency  dependent 
speckle  image  power  spectrum,  ^(F)  , in  accordance  with  the  exact 
definition  provided  by  Eq.  (18),  without  introducing  any  approximations 
relating  to  the  value  of  |'r-'r'|/r0  , or  to  the  value  of  X[F|/rQ  , in 

this  appendix  we  shall  carry  out  a numerical  evaluation  of  the  integral 
in  Eq.  (18)  for  various  values  of  D/rQ  and  of  \|f  |/rQ  using  Monte 
Carlo  methods.  This  allows  us  to  test  the  accuracy  of  the  results  we 
developed  in  Section  2 of  this  report  using  approximations  concerning 
the  values  of  |‘r-'r'|  / rQ  and  of  X 1 1* | / rQ  . The  numerical  results  we 
shall  develop  here  will  show  that  these  approximations  lead  to  quite 
accurate  results. 

Making  use  of  Eq.  's  (18),  (21),  and  (22),  we  may  write  the 
normalized  spectrum,  which  represents  the  speckle  transfer  function,  as 


J10)  = IT  W(r+xr)  W (r' ) W(f'+\T) 

X exp  {-6.88  r0-6'3[  |r-?'|B/&  + |xT|b/3 

- * (|r-?'  + Xf|R*  + . (A>1) 


It  will  prove  convenient  for  us  to  make  a change  of  the  variables  of 
gration  here.  For  this  purpose,  we  define  the  variables  x and  p 
ing  to  the  equations 

x = f/D 


inte- 

accord- 


(A.2) 
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from  which  it  follows  that 


dx  = dr/D» 

and 

?=(?-?')/D 

from  which  (treating  r asa  constant)  it  follows  that 
dp  = cT?'/D> 

We  can  now  rewrite  the  speckle  transfer  function  of  Eq.  (A.  1)  as 

j(o)  ~ J*J*  dx  dp  y(x)y(x+a)y(x-p)y(x-p  + a) 

x exp  f-6.  88(D/r0)B/3  [ |p|5*  + |£|a* 


- i (|  P + a|8/3 

+ |p-o|8*]} 

where 

y(x)  = < 

1"  1 » if 

1 

1*1  < i 

L o , if 

1*1  > * 

and 

a = \t/D 

* 

represents  the  normalized  spatial  frequency.  (The  effective  range  of  the 
normalized  spatial  frequency  is  from  zero  to  unity.  ) 

It  is  convenient  to  rewrite  Eq.  (A.  4)  as 

J(t) 

7(0)  ' JJdSdjSy®  VR  + S,  yR-jJ)  y(S.J+3) 

x exp  [-6.  88  (D/r0)B/3  ¥ (£,-£)] 
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(A.  2') 

(A.  3) 

(A.  3') 

(A.  4) 

(A.  5) 

(A.  6) 

(A.  7) 


where 


Y(p.o)  = |p|5/s  + |a|5/&  - | (|P  + aj6'5  + |p-a|B/3) 


(A.  8) 


Making  use  of  Eq.  (24),  we  note  that  accordingly  as  | a|  is  much 
greater  than  or  much  smaller  than  |p|  , we  can  approximate  Y(p,7)  by 

the  expressions 

l(p.a>=  |a|*«  - I |a!>|P|->»  f1  - 


if 


a « P 


(A.  9) 


and 


fp.3)  - |?r  - 1 !?l‘  |3|-l/3  1 - 


if  I a|  » |P| 


(A.  10) 


We  shall  not  use  these  approximations  in  an  analytic  way  to  develop  our 
results  (as  was  done  in  Section  2 of  this  report),  but  rather  will  use 
Eq.  's  (A.  9)  and  (A.  10)  to  provide  guidance  in  the  selection  of  the  sampling 
function  we  shall  use  for  the  Monte  Carlo  evaluation.  We  shall  also  use 
these  equations  in  our  numerical  computations  to  assure  adequate  accuracy 
in  the  evaluation  of  y(p,o)  when  the  magnitudes  of  |p|  and  | &|  are 
greatly  different,  without  resorting  to  double  precision  computer  compu- 
tation. 

For  Monte  Carlo  evaluation  of  the  integral  in  Eq.  (A.  7),  we  shall 
try  to  use  sampling  distributions  for  x and  p that  will  fully  sample  the 
range  of  allowed  values  of  x and  p , but  in  such  a way  as  to  minimize 
the  "variability  of  the  integrand.  " This  minimization  of  the  variability 
of  the  integrand  is  achieved  by  making  use  of  the  fact  that  for  Monte  Carlo 
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integral  evaluation,  we  can  multiply  and  divide  the  integrand  by  any 
desired  probability  density  and  then  "extract"  the  probability  density 
function  in  the  numerator  for  use  as  the  random  sampling  density  gen- 
erator. The  original  integrand  divided  by  the  probability  density  is  left 
as  a modified  integrand,  whose  mean  value  over  the  sampling  density  we 
wish  to  determine.  This  mean  value,  which  we  estimate  as  a Monte 
Carlo  average,  represents  the  value  of  the  integral.  The  accuracy  with 
which  the  Monte  Carlo  procedure  will  estimate  the  mean  value  of  this 
modified  integrand  depends  on  the  variability  of  the  integrand  in  the 
sampling  range.  If  we  hhve  been  reasonably  clever  in  our  choice  of 
the  sampling  distribution,  the  modified  integrand  will  have  a limited 
variability  and  it  will  take  relatively  few  samples  to  achieve  the  desired 
accuracy.  If  we  have  not  been  so  clever,  more  samples  will  be  required. 
Either  way,  the  validity  of  the  result  is  unaffected  — it  is  just  a matter 
of  how  many  samples  we  need  to  use  to  achieve  the  desired  accuracy  with 
our  Monte  Carlo  integral  evaluation  procedure.  In  general,  and  for  the 
problem  we  are  treating  here,  it  is  not  clear  what  practical  sampling 
distribution  is  optimum,  and  the  matter  of  making  an  efficient  choice 
is  somewhat  of  an  art.  We  can  judge  our  success  in  making  a suitable 
choice  of  the  sampling  probability  density  by  evaluating  not  only  the 
mean  value  but  also  the  standard  deviation  of  the  modified  integrand 
when  we  generate  our  Monte  Carlo  average,  and  in  fact  we  can  use  this 
standard  deviation  to  help  us  decide  when  we  have  taken  enough  random 
samples. 

For  Monte  Carlo  evaluation  of  the  integral  in  Eq.  (A.  7),  we  chose 

to  use  a uniform  distribution  for  x and  a gaussian  distribution  for  p . 

We  chose  the  range  of  the  x- distribution  to  correspond  to  the  rectangle 

that  just  encloses  the  region  of  overlap  of  the  two  circles  defined  by 

y(x)  and  y(x  + a)  . (Outside  this  rectangle,  the  integrand  has  zero  value 

* It  will  later  become  clear  from  a study  of  Table  A.  2 that  we  have  not  been 
very  clever  in  our  choice  of  the  sampling  distribution  when  D/rQ  is 
very  large  and  lal  **rQ/D  • 
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so  nothing  is  excluded  by  use  of  the  limited  sampling  range. ) If  we  denote 
the  component  of  x that  is  parallel  to  j by  Xj  , and  the  perpendicular 
component  by  Xg  , then  our  rectangle  corresponds  to  the  region  defined 
by  the  equations 

-fc  < < £ - | »|  • (A.  11a) 


where 


-ih  **  < h] 


h = (1  - |2|8FS 


(A.  lib) 


(A.  12) 


This  is  a rectangle  of  width  (1  - |al)  and  of  height  (h)  . The  probability 
density  function  is  thus 

[ (1-  iJlHO-i.  i{  (-4<*i<jl«l) 

and  {-^h<Xg<^h}  (A.  13) 

k 0 , otherwise 


As  remarked  above,  we  will  use  a (two-dimensional)  gaussian 
probability  density  foT  the  p sampling  distribution.  We  can  write  this 


P„  (?)  = (2tt  exp  (-£  Jp| */oa) 


(A.  14) 


where  a is  a sampling  parameter  whose  value  we  have  to  choose.  To  get 
some  insight  into  how  to  choose  a value  for  a , we  multiply  and  divide 
the  integrand  of  Eq.  (A.  7)  by  Py  and  by  Ps  . Thus  we  have  for  the 
speckle  transfer  function 

= ($")■*  fj  dx  d?y(x)  y(x+ a)  y(x-p)  V(x-p  + a) 


P0tf)  P6<P) 

X exp  [-6. 88  (D/r0)B^  y(p,a)]  p • 


(A.  15) 
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Making  use  of  Eq.  's  (A.  13)  and  (A.  14),  this  can  be  rewritten  as 

~ °*  h(l-|3|)  JJd3<ipP„(S)  Pc(p-) 

X {>£)>£+ a)  *£-?)  >£-?+£) 

X exp  [*  |p|»/a*  - 6.88  (D/rQ)s*  Y(p.S)]}  . (A.  16) 

The  quantity  in  the  large  curly  brackets  is  the  modified  integrand  whose 
mean  value,  in  accordance  with  the  probability  densities  Pu  (x)  and  P (j$), 
we  wish  to  determine  by  Monte  Carlo  random  sampling. 

Taking  note  of  Eq.  (A.  10),  we  reason  that  we  can  minimize  the 
effective  variability  of  the  modified  integrand  when  Ja|  is  mudi  larger 
than  r0/D  by  choosing  the  sampling  parameter  a to  make  the  argument 

of  the  exponential  in  the  modified  integrand  vanish  when  |p|  equals  0 

This  leads  to  the  condition 

a = (r0/D)  (£/6.  88)^ 

0.207  (r0/D)  , for  |a|  » (rQ / D)  . (A.  17') 

If  we  now  take  note  of  Eq.  (A.  9),  we  see  that  when  |£|  is  much  less  than 
rQ/D  , then  the  argument  of  the  exponential  in  the  modified  integrand  will 
have  relatively  little  variation  with  p if  c is  comparatively  large.  So 
as  to  avoid  taking  too  many  samples  outside  the  range  in  which  y(x-p)  and 
^(x-p  + a)  are  non- zero,  we  somewhat  arbitrarily  suggest  that  the  sampling 
parameter  ^ have  the  value 


a = * ’ for  I « (r0/E»  • (A.  17~) 


To  smoothly  combine  Eq.  's  (A.  17')  and  (A.  17"),  we  chose  to  use  as  the 
gaussian  distribution  sampling  parameter 


I 


(A.  17) 


Given  a , it  is  a straightforward  matter  to  generate  appropriately 
distributed  random  samples  of  the  two  components  of  p , namely,  px  and 
ps  according  to  the  equation 


Px  - a gx  ; P2  = a g2 


(A.  18) 


where  gl  and  g3  are  gaussian  distributed  random  variables  with  zero 
mean  value  and  unity  standard  deviation.  Well-defined  computer  routines 
for  generating  the  random  variables  gt  and  gg  exist.  The  generation 
of  the  random  samples  Xj  and  Xg  follows  from  the  equations 

*i  - |a|  - (!-|a|)  Ul  , (A.  19a) 


S = h (u2  - ^) 


(A.  19b) 


where  Uj  and  ug  are  random  variables  uniformly  distributed  on  the 
interval  between  zero  and  unity  (0,1).  [Clearly  xJ  and  Xg  are  uniformly 
distributed,  and  by  substituting  the  limiting  values  of  u and  u2  , we  can 
confirm  that  xt  and  just  span  the  ranges  defined  by  Eq.  's  (A.  11a)  and 
(A.  lib).  ] 

If  we  use  the  notation  MCA  {jf x;  pt  (x)  | y;  pg  (y)  ||  f(x,y)J  to  denote 
the  average  value  of  the  function  f(x,^)  when  evaluated  as  a Monte  Carlo 
Average  over  the  sampling  distributions  of  x and  y defined  by  px  (x)  and 
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Pa  (y)  » then  it  is  obvious  from  consideration  of  Eq.  (A.  16)  that  we  can 
write  the  speckle  transfer  function  as 

™ a»h(l-|3|)MCA{*';Pu<2)|?;P#6)| 

(y(x)  > (x  + a)  >(x-p)  >(x-p  + 

X exp  [fc  |?|  /a*  - 6.  88  (D/rQ)e* 

We  can  estimate  the  variance  of  the  modified  integrand  as 
Var  ~ ” a8  h (1-|3|)  MCA  {3 ; P„  (2)  |? ; PQ  (?)  || 

{y(S)  >£+£)  ytf-p)  ytf-?+£) 

X exp  [fc  !p|  /o3  -6.  88  (D/rQ)s*  *(p,2)]38  j 

- [j(T)//(0)3> 

Accordingly,  the  standard  deviation  of  the  Monte  Carlo  estimate  of  the 
normalized  bispectrum,  j(£)lj(0)  , which  we  obtain  from  Eq.  (A.  20) 
is 

SD  = (Var/C)i/a  , (A.  22) 

where  C is  the  number  of  sample  values  of  (x,?)  that  we  used  in  forming 
the  Monte  Carlo  Average. 

In  Appendix  B,  we  list  a computer  program  which  generates  a table 
of  values  for  the  speckle  transfer  function,  0)  , and  for  the  associated 

standard  deviation,  subject  to  the  conditions  that  no  more  than  one  million 
samples  will  be  used,  and  that  otherwise  the  Monte  Carlo  sampling  process 


l 


(A.  21) 


(A.  20) 


r 


y- 


a* 


I 


will  stop  as  soon  as  an  estimated  one -percent  accuracy  for  the  speckle 
transfer  function  is  achieved.  Results  were  calculated  for  values  of 
D/rQ  equal  to  10  , 30  , 100  , 300  , and  1000,  and  for  a very  large  range 
of  values  of  the  normalized  spatial  frequency,  a = \T/D  , from  very  near 
zero  to  very  near  unity.  The  choice  of  valuej  for  the  magnitude  of  a was 
made  to  match  a gaussian  distribution  scale  so  as  to  provide  fine  detail  for 
|o|  near  zero,  as  well  as  for  |al  near  unity.  The  speckle  transfer  func- 
tion values,  0)  , are  listed  in  Table  A.  1 . In  Table  A.  2,  we  show  the 

ratio  of  the  standard  deviation  to  the  estimated  value  along  with  the  number 
of  samples  used  to  form  that  estimated  value  for  each  entry  in  Table  A.  1.  In 
evaluating  the  speckle  transfer  function,  random  samples  were  used  in  blocks 
of  one-thousand,  the  number  of  blocks  used  being  adjusted  to  provide  one- 
percent  accuracy,  except  that  no  more  than  one  thousand  blocks  were  ever 
used. 

We  recall  that  our  objective  in  this  appendix  was  to  provide  a basis 
for  testing  the  accuracy  of  the  approximations  to  the  speckle  transfer  func- 
tion provided  by  Eq.  (40)  for  large  values  of  |a|  relative  to  rQ/D  .and 
by  Eq.  (44)  for  small  values  of  |o|  relative  to  rQ/D  . Accordingly,  our 
computer  program  also  calculated  these  approximation  results.  In  Fig.  's  1 
to  5,  we  have  plotted  the  Monte  Carlo  estimates  of  the  speckle  transfer 
function  as  a function  of  the  normalized  spatial  frequency,  | = X |T|/D  » 

for  each  of  the  five  values  of  D/rQ  considered,  along  with  the  correspond- 
ing approximation  results  of  Eq.  's  (40)  and  (44).  As  can  be  seen  from  a 
consideration  of  these  figures,  the  approximations  are  each  quite  good  in 
their  region  of  expected  validity.  This  is  particularly  so  for  the  larger 
values  of  D/rQ  . 
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Table  A.  1 


Speckle  Transfer  Function 
Monte  Carlo  Results 


Normalised 


Spatial 


SpacUa  Tra 


afar  Function 


Frequency 

lal 

D/r0  • 10 

D/r#  • 30 

D/r0  - 100 

D/ r0  « 3 00 

D/r.«  1000 

C. 00021 

9.9t«-0l 

9.92, 

-Cl 

9.81, 

01 

9.29# 

-Cl 

4.14,-Cl 

z.zzozt 

i.:s«- Jo 

l.CI# 

-CO 

9.67# 

01 

8.87, 

-01 

3.01,-01 

0.05012 

l.CW-30 

9.93, 

-01 

9.81# 

01 

9.59# 

-01 

1.72.-01 

c. 030*0 

9.98,-31 

9.95, 

-01 

9.65# 

01 

8.27# 

-01 

2.44#-fll 

c.'coeso 

9.ee*-:i 

9.89, 

-01 

9.56# 

01 

7.62# 

-01 

1 .14,-Cl 

3.020*2 

9.89,-31 

9*  9 3t 

-01 

9.14, 

31 

6.76# 

-01 

5.71,-02 

C .02077 

9.80,-31 

9 . 9 1 • 

-01 

9.C?, 

01 

5 . ?!» 

-01 

i 7*,-or 

C . COO  96 

».96t-3i 

9.81, 

-Cl 

9.83# 

01 

4.62# 

-01 

1.15, -Cl 

o.r?: i9 

9.96,-31 

9.72* 

-c: 

8.25# 

01 

1.2', 

-c: 

2.80,-04 

C.C3147 

9.96,-31 

9.63# 

-Cl 

7.57, 

01 

2.06, 

-Cl 

1. 18,-CS 

2.021*2 

9. 78,-31 

9.42* 

-Cl 

6.88# 

51 

1.13, 

-Cl 

1.21.-C6 

:.co?i5 

9.«4*o: 

9.24. 

-Cl 

5.95# 

31 

4.47# 

-It 

7.54,-07 

1.02279 

9. 51,-01 

8.95, 

-01 

4 • 9 C* 

31 

1 .26, 

-Cl 

8.95,-07 

9.59#-:: 

8.56# 

-o: 

1.61, 

ji 

2.32, 

- y 

7.84,-07 

-.rctii 

9.52*-:i 

9.25# 

-Cl 

2.50, 

31 

2.40# 

-0* 

7.61,-07 

c,r?f  is 

5.54*-31 

7.63# 

-o: 

1.48# 

31 

2.69# 

-c? 

6. 80#-07 

9. ?3*-31 

6.87# 

-01 

7.  12# 

32 

1.28. 

•OS 

6.91,-07 

0.2:77* 

9.52,-01 

5.99# 

-01 

2.71* 

J2 

1 .02. 

-05 

6.43,-07 

o.c:**3 

1.74,-31 

4.92# 

-Cl 

7.61, 

33 

9.25# 

-C6 

*.05#-27 

o.c: :*4 

S.SH-31 

3.90# 

-01 

1.58# 

33 

e.50# 

-C6 

6. 06#-07 

c.c: ?sw 

*•019-31 

2.89# 

-Cl 

1.10, 

34 

C.27# 

-C6 

5. 73,-07 

0.0!ta6 

76#- 3 1 

1.91# 

-c: 

1.32, 

34 

7.58. 

-06 

5.67,-07 

0.02C32 

6.  62#- 31 

1.17# 

-01 

1.04# 

04 

7.47, 

-0* 

«. 35,-07 

C.02?9* 

5.85,-51 

5.75# 

-02 

8.89# 

05 

7.03. 

-C6 

5.*2,-C7 

0.02649 

5. 11, -01 

2.61# 

-C2 

8.22, 

25 

6.91# 

-06 

5.20,-07 

0.P1177 

4.26#-31 

1.05# 

-02 

7.61, 

05 

6.50* 

-0. 

5.22,-37 

G.C39!? 

1.11,-31 

4.21# 

-03 

7.17# 

05 

6.46# 

-06 

5.31.-C7 

0.0469? 

2 • SO*- 31 

2.10, 

-03 

6.84* 

35 

6.15* 

-06 

4.95.-07 

0*05464 

i.9o#-oi 

1.19* 

-03 

6.78, 

33 

6.01* 

-Cb 

4. 9*. -07 

0.06353 

1 . 7 0# - 3 1 

1.10# 

-03 

6.25, 

05 

5. 90* 

-06 

4.72*-C7 

0. 07550 

8.42,-32 

9.66# 

-04 

5.91, 

05 

5.77# 

-06 

4.88,-07 

O.094M 

5.3b#-32 

8.(1, 

— C 4 

5.79# 

55 

5.53. 

-06 

4. «?»-D? 

0.09691 

7 • 50«- 32 

8.09# 

-04 

5.51# 

35 

5*43# 

-06 

4. 36#-27 

o. 1! 046 

2.24#-02 

7.57# 

-04 

5.35# 

35 

5.26# 

-06 

4 • 23*-C7 

0.12529 

1*6  0 # - 32 

7.15* 

-C  4 

5.0?, 

25 

4.99# 

-06 

4.19,-07 

0.14142 

1. 20# -32 

6.69# 

-04 

4 • 86# 

05 

4.86# 

-06 

1.59,-07 

o.  15888 

9.9*«-03 

6. 3-# 

-04 

4.67# 

05 

4.63* 

-G6 

1.81,-07 

0.17766 

?. 33#-33 

5.95# 

-04 

4.43# 

35 

4.45* 

-Cb 

1.74#-C7 

0.1*777 

7. 38#-33 

5.67# 

-04 

4.20# 

05 

4.27* 

-06 

?• 65#-G7 

0.21919 

: • t5*-33 

S.42# 

-C4 

4.08# 

35 

4.12# 

-06 

1.48,-37 

o.2-i a? 

7.C 3#-33 

5.09# 

-04 

3.78# 

05 

3.65# 

-06 

I. 11,-07 

C.  26561 

c . b-*-33 

4.83# 

-04 

3.66# 

35 

3.70# 

*06 

I. 19#-27 

0.2?0?4 

5. 39#-: J 

4.43# 

-0a 

3.46# 

35 

3.44* 

-06 

3 • D2#-C7 

0.71716 

4 • 5 3*-3  3 

4.23# 

-04 

7.20# 

35 

3.30* 

-06 

2.87,-07 

0.3897C 

? 3 

7.54# 

- C 4 

2.839 

05 

2.95* 

-06 

2.51,-07 

0.42265 

3. 77*- 33 

3. 06* 

-0- 

2.50* 

35 

2.56* 

-06 

2.26,-07 

0.47826 

2.7H-:j 

2.62# 

-c* 

2.11, 

05 

2-20. 

-06 

1.94#-C7 

3.53764 

?.25#-:j 

2.16# 

-04 

1.30# 

35 

t .94* 

-06 

1.63*-27 

3.53373 

l. 79*-;3 

1.77# 

-24 

1.50# 

35 

1.57* 

-06 

t . ?9#-G  7 

3.64263 

1.45.- J3 

1.42# 

-24 

1.23, 

33 

1.25. 

-06 

1.11,-07 

0.69145 

t .09fc-:3 

1 • l 5* 

-04 

9.81# 

36 

1.01* 

-06 

9 • 94  #-0  8 

0.74067 

3 • 37t - : 4 

8.97, 

-05 

7.55, 

06 

8.03. 

-07 

6.92#-C8 

C. 73574 

5 • 16#- j4 

6.79# 

-25 

5.64# 

36 

6.23* 

-07 

5.11,-08 

C.  822  54 

4. 36t- 

4.9  7# 

-C  5 

4.31# 

06 

4.54* 

-07 

4.04*-Cf 

0.85613 

1.01,-;. 

3.60* 

-03 

1.10* 

36 

1.35. 

-07 

2.02*-:! 

0.63571 

2. 02*- 34 

2.44# 

-05 

2.21# 

06 

2.39. 

-07 

2 . 08#-39 

0.91062 

1 • 23«- 34 

1.64# 

-05 

1.53# 

36 

1.66* 

-07 

1 .44*-C3 

0.93129 

7. 40«-3  5 

1.07# 

-05 

1.  JO, 

36 

1.106 

-07 

9. 77#-39 

0.94839 

4. 01«- 35 

6 • 8 C# 

-Ce 

6.51# 

37 

7.34* 

-C9 

*.5J*-:9 

0.961-5 

2 , 0 1 #- 3 5 

7.93* 

-C6 

4.U# 

27 

4 . 5**# 

-29 

4.026-39 

0.97133 

9.46,-36 

2.21# 

-38 

2.47# 

37 

2.87, 

-08 

2.6!«-C9 

0.97974 

4.18,-06 

1.119 

-26 

1.45# 

37 

1.69* 

-OS 

1 .55*-C9 

0.98563 

1.75#- 36 

*.34* 

-37 

5.39* 

29 

1 .01# 

-03 

o.26*-lC 

0.99994 

8.S8#-37 

:.3»« 

-C  7 

4.75# 

39 

5.77. 

-09 

5.27,-10 
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Monte  Carlo  Error  Parameters 

For  each  of  the  Monte  Carlo  results  shown  in  Table  A.  1 we  show  here  our 
estimate  of  the  rms  fractional  error,  a , and  the  number  thousands  of 
samples,  C , required  to  achieve  this  fractional  error.  Because  we  never 
allowed  use  of  more  than  one-million  samples  (C=1000),  our  worst  case 
fractional  error  was  19.  5%  (SD=0.  195),  but  for  the  most  part  a 1%  (SD=0.  01) 
fractional  error  was  achieved,  generally  with  only  about  37,000  (C=37) 
samples  in  the  Monte  Carlo  estimation  process. 


Normalised 

Spatial 

Frequency 

D/r#- 

10 

D/r#« 

10 

D/r#. 

100 

D/r0  = 

100 

D/r,  * 

1000 

1*1 

SO 

c 

SO 

c 

SD 

c 

SO 

c 

SD 

c 

0*00021 

S.010 

17. 

0.  313 

C.310 

17. 

0.010 

11. 

9 . C 1 9 

17 

0.0CC26 

C.310 

17. 

o.:*.: 

16. 

0.310 

57. 

o.cio 

11. 

9.  Cl  9 

37 

0.00C12 

0.310 

17. 

0.910 

!•  . 

O.CIC 

J7. 

0.010 

57. 

9.910 

5* 

0.00040 

0.010 

16. 

0.31C 

p. 

3.310 

37. 

0.010 

17. 

0.313 

16 

0*00050 

C.01C 

i*. 

o.oi: 

J 7 « 

0.010 

57. 

C.010 

37. 

3.310 

36 

0.03CS2 

0.310 

i*. 

0.013 

37. 

:.cia 

17. 

O.OIO 

36. 

onto 

55 

0*00077 

5.313 

17. 

o.cic 

17. 

:.3io 

17. 

0.010 

37. 

9.013 

55 

0*00096 

0.010 

17. 

0.01C 

37. 

3.  010 

17. 

o.oio 

16. 

3.910 

55 

0.00119 

0.9*.  9 

17. 

0.010 

I*. 

0.013 

17. 

0.010 

36. 

?.  313 

* 5 

0.03147 

9.310 

17. 

0.010 

17. 

3.  310 

16. 

0.010 

16. 

0.  044 

1090 

c.oom 

0.310 

17. 

o.oio 

17. 

3. CIO 

17. 

3.010 

1*. 

0.  1 65 

1339 

0.00225 

0.313 

56. 

o.oic 

i’. 

3.  313 

17. 

3.310 

32. 

3.  147 

1000 

3.03276 

3.010 

!6. 

O.CIC 

!?. 

3. CIO 

17. 

0.010 

31. 

0.  117 

1900 

0. 003-3 

3.313 

17. 

0.010 

• : . 

3.  010 

16. 

0.010 

47. 

3.  ?9 4 

1930 

0.03422 

3.31C 

17. 

Q.31G 

M. 

:.on 

15. 

0 . C 1 3 

911. 

0.  065 

IS  00 

0.0CSI* 

3.310 

17. 

0.313 

3.913 

54. 

3.054 

1 9C9. 

0.244 

1090 

0 .006  34 

3.310 

37 . 

x - . 

5 3. 

0.066 

1 COS. 

0.011 

I 33  2 

0.03774 

0.01  0 

17. 

0.010 

J 7 • 

:i: 

32. 

3.C52 

1 503. 

3.021 

10  30 

C. 00941 

0.010 

17. 

o.si: 

36. 

:.  :io 

16. 

3.040 

1333. 

3.015 

1930 

0.011U 

0.310 

17. 

3.  CIO 

54. 

0.01  3 

235. 

0.C29 

1303. 

0.313 

1C3Q 

0.01334 

:.3ic 

38. 

0.013 

34. 

:i5 

:::c. 

3.323 

1999. 

3.  313 

464 

0.01666 

C.3IC 

u. 

0.310 

3 J. 

:.024 

13  30  . 

C.314 

1 990. 

3 . 0 ! 3 

215 

0.02002 

0.310 

i». 

0.310 

ii. 

j.  321 

1030. 

3.310 

1 903. 

0.012 

114 

0. 02194 

5.910 

17. 

0.3 1C 

12. 

:i7 

1 C 30. 

9.019 

515  • 

9.019 

53 

C.  02949 

9.010 

I*. 

0.010 

61. 

3.::i 

1030. 

C.010 

270. 

9.010 

27 

9.05577 

0.910 

1*. 

3.01C 

S3. 

o.  no 

51*. 

o.cic 

143. 

0.013 

11 

C. 03993 

0.010 

1*. 

o.oic 

236. 

3.313 

514. 

3.910 

71. 

9.913 

6 

0.04677 

0.310 

16. 

9.010 

521. 

3.913 

2 59. 

0.010 

19. 

3.009 

4 

0.g;*»« 

3.010 

34. 

0.010 

62  S. 

3.013 

157. 

3. 010 

20. 

3.306 

4 

0.06353 

0.  010 

34. 

0.310 

.?». 

0.310 

?0. 

3.310 

13. 

9.99? 

9 

0. 07553 

9.310 

37. 

0.010 

141. 

3.  319 

52. 

3.0  09 

6. 

3.334 
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